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FOREWORD 

This report presents the results of a study performed by 
personnel in the Fluid Technology Application Section of the 
Lockheed-Huntsville Research & Engineering Center in fulfill- 
ment of the requirements of Contract NAS8-33801 with NASA- 
Marshall Space Flight Center. The NASA-MSFC contracting 
officer's tecluiical representative for this study was Joseph L. 
Sims, ED3 3. 
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Symbol 


A 

V A * 



Greek 

y 

6 


v 

0 


SYMBOLS AND NOTATION 


Description 


. 2 

area, in 

nozzle area ratio, dimensionless 


diameter, in. 

Mach number, dimensionless 

nozzle exit plane Mach number, dimensionless 

mass flow rate, Ibm/sec 

2 

static pressure, Ibf/in 

2 

vehicle base region pressure, lbf/in 

2 

nozzle chamber pressure, lbf/in 
nozzle exit plane static pressure, lbf/in^ 
jet boundary static pressure, lbf/in* 1 
nozzle chamber temperature, R 


ratio of specific heats, dimensionless 

initial plume expansion angle measured from the nozzle 
centerline, deg 

Prandtl- Meyer turning angle, deg 
nozzle wall angle, deg 


Subscr 

b vehicle base region conditions 

c nozzle chamber conditions 

e nozzle exit plane conditions 

j jet boundary conditions 

oo wind tunnel freestream conditions 

Supcrscr ipt 

❖ nozzle throat conditions 
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1. INTRODUCTION AND SUMMARY 

This study was conducted in support of an experimental program to 
define power -on aerodynamic characteristics of the Space Shuttle over a 
range of ascent trajectory conditions. 

MSFC has developed procedures to simulate Space Shuttle launch vehicle 
main propulsion jet plumes for use in aerodynamic wind tunnel tests using air 
to simulate the rocket plume. These procedures will be used by Rockwell 
International/Space Systems Group (RI/SSG) in conducting power-on aerody- 
namic wind tunnel tests of the thrust augmented Space Shuttle launch vehicle. 
These procedures include similarity parameter definition, data analysis pro- 
cedures and model nozzle design criteria to match values of the prototype 
similarity parameters with the auxiliary air supply available at the wind tunnel. 

The present form of the similarity parameter for base pressure match- 
ing is 

M. 6. \ 

. J . i. \ 

. . a b I 
M y. / 

e J/ Prototype 


M. 6. 

. . J. J 
a b 

M y . , 

c J / Model 


where the values of a and b are functions of the freestream Mach number. 
Mj is the Mach number on the plume boundary, 6 ^ is the initial expansion 
angle of the plume, y. is the ratio of specific heats of the rocket exhaust gas 
at the plume boundary and M g is the Mach number of the flow in the nozzle 
exit plane at the nozzle lip. This form of the similarity parameter was 
developed by the correlation of experimental base pressure data from a 
number of model nozzles flowing Freon 14 (C.F 4 ) as a prototype gas and air 
as the simulant gas. No implicit similarity conditions were used in this 
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correlation. This means that no single physical characteristic of the proto- 
type engine needs to be matched by the model nozzle. Thus, for example, 
it is permissible to change the nozzle exit wall angle in order to change 6j 
to achieve the correct numerical value of the similarity parameter with the 
model nozzle . 

Sections 2 and 3 of this report present the computational steps of an 
analytical procedure used to design an "envelope" of model nozzles (the 
envelope is defined by nozzle exit plane Mach number, M , and nozzle exit 
wall angle, l> e ) which may be used to simulate Space Shuttle launch vehicle 
main propulsion jet plumes over a specified range of ascent trajectory condi- 
tions. Each nozzle in this "envelope" will be capable of achieving the desired 
similarity parameter values within the restriction of a fixed maximum auxiliary 
air pressure supply and mass flow rate imposed by the facility where the wind 
tunnel tests will be conducted. The analytical procedure, though developed in 
support of a Space Shuttle experimental program, is not limited to Space Shuttle 
vehicle configuration applications. In fact, the analytical procedure may be 
used to design model nozzles to simulate the prototype engine conditions of 
any single body-single nozzle; single body-triple nozzle or triple body engine 
system (see Fig. 3.3 for a definition of these engine systems). 

A description of a computer code which automates the analytical pro- 
cedures used to design the above "envelope" of model nozzles is presented 
in Section 4 of this report. An example of the computer code input require- 
ments and resultant outputs are presented for the wind tunnel test IA-604 
planned for the NASA- Ames Research Center's 11x11 wind tunnel (Ref. 1-1). 
Final model nozzle designs and power sweep operating characteristics for 
each of the desired flight conditions to be simulated in the wind tunnel tests 
is presented for representative model nozzles that exist in the "envelope" 
of designed model nozzles. 
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2. OUTLINE OF THE ANALYTICAL PROCEDURE 


Tins section presents an outluu* of the computational steps of a com- 
puter code uned to design mode! n. ^zies which may be used to simulate the 
prototype engine conditions of any sing!. .* body -single nozzle, single body-triple 
triple nozzle or triple body engine system. Eacn model nozzle will be cap- 
able of achieving the desired similarity parameter values within the restric- 
tions of maximum air pressure supply and mass flow rate imposed by the 
facility where wind tunnel testing will be performed. Section 3 of this report 
presents a more detailed discussion of each computational step performed 
by the computer code as well as a discussion of the computer code input re- 
quirements. It is strongly suggested that the computer code user make every 
possible use of Sections 2 and 3 of this report when setting up his data deck. 
Doing so will give the user a greater understanding of the model design 
problem and the actual use of the computer code. 

The basic equations employed by the computer code are given in Table 

2 - 1 . Examination of these equations indicates that the design problem can 

be resolved by considering three primary variables to define the candidate 

model nozzle designs. These are: (1) model nozzle chamber pressure (P ); 

>i< c 

( 2 ) model nozzle exit plane Mach number (M ) or area ratio (Ae/A ); and 

(3) model nozzle exit plane wall angle (0^). 

Utilizing the nozzle flow analysis described in Ref. 2-1 for a gaseous 
only flow, or Ref. 2-2 for a gas-particle flow and an initial plume Prandtl- 
Meyer expansion calculation, a parametric variation of the above three 
parameters is performed to define candidate model nozzle solutions which 
satisfy the base pressure similarity parameter relationship 



model 


- M. 6. ' 

— 1 - 4 - 

M a r b 
* e J 


prototype 
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T a b U‘ i-t 

BASIC EOLATIONS 1'SED IN THE: MODEL 
NOZ/LK DESIGN ANALYSIS 


The nozzle design analysis is performed assuming lhat the model 
no/zli will be flowing a constant ) simulant gat during the wind tunnel 
testa, therefore; 




mode I 


> } > constant 

c e J 


/M.£ \ 
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' l> I ' 
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The final result is an envelope of candidate model nozzles (the envelope is 
defined by M > and G^) which may be used to simulate the flight conditions 
that exist for all flight Mach numbers to be tested. The wall contour of the 
actual model nozzle used in the wind tunnel test is arbitrary and left up to 
the nozzle designer. Within the scope of the current technology it is felt 
that a conical nozzle with a. designed value of M^ and 6^ will suffice for all 
cases . 


• COMPUTATIONAL STEPS 

The specific steps of the analytical procedure used to design model 
nozzles are outlined below. With the exception of actually calculating the 
value of similarity parameter in step 4, the first 6 steps of the analysis are 
performed by the computer code user and input to the computer code in the 
form of a data deck (see Section 4.2 for a description of the computer code 
input requirements.) The remaining computational steps have been auto- 
mated by the computer code. 


1 . A schedule of flight Mach numbers that is to be simulated is 
determined as a function of the test purpose. 

2. For each flight Mach number to be simulated, the correspond- 
ing values of freestream static pressure, Per* and prototype 
nozzle chamber pressure, P , are determined. 

i. A predicted value of base pressure ratio (Pb/P^) 8et f° r 
each flight Mach number to be simulated. The value of P b/ Poo 
used is based upon experience gained from previous wind 
tunnel tests . 

4. Prototype plume similarity data are determined. These data, 
together with the above P^/ P^ and P c , are used to determine 
the value of prototype base pressure similarity parameter to 
be simulated to ensure that base pressure matching occurs. 

5. The physical restrictions of prototype nozzle exit diameter, 
scale of the model nozzle, maximum pressure supply and 
mass flow rate imposed by the facility are set. 

6. A gasdynamic constraint on the ratio P e /P^, is set to ensure 
that model nozzle flow separation does not occur. The minimum 
allowable value of P ^P^ determined experimentally is 0.60. 
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7. Scliil a flight Mach number from the schedule of flight Mach 
mnuhci’tt to be simulated . 

8. Select a value of model nozzle exit plane Mach number. The 
value of M t , will be varied parametrically in subsequent itera- 
tions. Selecting a specific val^e of Mach number gives the 
model nozzle area ratio, A /A . 

9. The next design parameter of concern is model rvazzle chamber 
pressure, P. , Selecting a specific value of P„ then permits the 
calculation of jet or plume boundary Mach number, Mj . The 
value of P c chosen is usually set equal to 50% of the maximum 
available air pressure supply. 

10. The value of 6j can be found which ensures similarity using the 
following relationship: 


/model 



model 



prototype 


11. Finally, (he model nozzle exit plane wall angle, 0 , required to 

achieve A A is then 

•'/model 


6 i) * *0 ‘ 

■V mode l /model '/model 

Incrementing the value of M and repeating steps 10 and 11 will 

define a family of candidate nozzles which will; (1) ensure base 
pressure matching for 'one" specified set of flight conditions and 
satisfy the mass flow rate and maximum air pressure supply 
constraints set by the test facility. 

1Z. Steps 7 through 11 are repeated for each of the flight Mach num- 
bers to be simulated. 

1 1. Once step 1Z is completed, each family of candidate model nozzles 
derived for a particular flight Mach number is checked to see which 
nozzles are capable of simulating the flight conditions that exist for 
all other flight Mach numbers. 
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14. The final result is an envelope of model -or.zles (the envelope is 
defined by and ©^) which may be used l. simulate the flight 

conditions that exist for all flight Mach numbers to be tested. 

15. Calculate the power sweep operating characteristics for each 

of the desired flight conditions to be simulated r or representation 
model nozzles that exist in the envelope. 


A flow chart of a computer code used to des gn model nozzles which 
meet MSFC base pressure similarity parameter criteria is presented in 
Chart 2-1 . 


2 -5 


lOCKHEFn HUNTSVILLE RESEARCH & ENGINEERING CENTER 















I 

s 

> 

I 


tt 









LMSC-HREC TR D7841 1 1 


3. A DETAILED DESCRIPTION CD THE COMPUTATIONAL STEPS 
OF A COMPUTER CODE USED TO DESIGN MODEL NOZZLES 
WHICH MEET MSFC RASE PRESSURE 
SIMILARITY PARAMETER CRITERIA 

This section presents a detailed discussion o£ each computational step 
as well as input requirements of a computer code used to design model nozzles 
used to simulate the prototype engine conditions that exist over a specified 
range of ascent trajectory conditions for any single body- single nozzle, single 
body-triple nozzle or triple body engine system. The general engine systems 
are defined in Fig. 3-3. Each model nozzle will be capable of achieving the 
desired similarity parameter values within the restriction of maximum air 
pressure supply and mass flow rate imposed by the facility where wind tunnel 
testing will be performed. 

It is strongly suggested that the computer code user make every possible 
use of Sections <1 and 3 of this report when setting up his data deck. The user 
will thus gain a greater understanding of the model nozzle design problem and 
the actual use of the computer code. 

3.1 FUNDAMENTAL EQUATIONS 

The basic equations employed by the computer code are presented 
on the following page. 

The nozzle design analysis is performed assuming that the model 
nozzle will be flowing a constant y simulant gas during the wind tunnel 
tests, therefore 

y model 1 >"c = V e = y , = con8tant < 3 • 1 > 
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Solving for M in Eqs.(3.2) and (3.3) yields 



(3.6) 


The equation for plume boundary initial angle is 


6 . = 

J 


0 + 1 / 


(3.7) 


and solving Eq. (3.7) for 0 e yields 


0 

e 


6 . + v - V- 
J e J 


(3.8) 


Base pressure matching is achieved when 


Similarity Parameter) , . = Similarity Parameter) t . 

/ model 7 J prototype 


or, 


M. 6 . 

J. J 

l M a y b 
k e j 


model 


M. 6. 

x . a ,b 
M y . 
e J 


(3.9) 


prototype 


(See Symbols and Notation for description of variables.) 
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Examination of Eqs.(3.1) through (3.9) indicates that the design problem 

can be resolved by considering three primary variables to define the candidate 

model nozzle designs. These are: (1) model nozzle chamber pressure (P ); 

% c 

(2) model nozzle exit plane Mach number (M ) or area ratio (A /A ); and (3) 
model nozzle exit plane wall ai.gli (O '. 

Utilizing a nozzle flow analysis described in Ref. 2-1 for a gaseous 
only flow, or Ret . 2-2 for a gas-particle flow and an initial plume Prandtl- 
Meyer expansion calculation, a parametric variation of the above three 
parameters is performed to define candidate model nozzle solutions which 
satisfy the base pressure similarity parameter relationship 


M. 6. 

-J— JL 
u a, b 
M } . 

1 j / model 



prototype 


The final result is an envelope of model nozzles (the envelope is defined by 

M^ and (< ) which may be used to simulate the flight conditions that exist for 

all flight Mach numbers to be tested. The wall contour of the actual model 

nozzle used in the wind tunnel test is arbitrary and left up to the nozzle 

designer. Within the scope of the current technology it is felt that a conical 

nozzle with a designed value of M and 0 will suffice for all cases. 

e e 

3.2 SIMILARITY PARAMETER DEFINITION AND APPLICATION 

Before proceeding further in the description of the analytical design 
techniques a better understanding of the similarity parameter and how it is 
used in the data analysis stage of the test program is necessary. 


The similarity parameters used in this analysis have been supplied by 
MSFC. The present similarity parameter (SP) used for base pressure match- 
ing is of the general form 


M. 6- 

SP) . - — J-4- 

gene ral M a b 

e j 


(3.10) 
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where the values of a and b are functions of frecstream Mach number and 
vehicle test configuration. This form of the similarity parameter was devel- 
oped by the correlation of experimental base pressure data from a number of 
model nozzles flowing Freon 14 (CF^) as a prototype gas and air as the simu- 
lant gas. No implicit similarity conditions were used in this correlation. 

This means that no single physical characteristic of the prototype engine 
needs to be matched by the model nozzle. Thus, for example, it is permiss- 

able to change the nozzle exit wall angle in order to change 6. to achieve the 

J 

correct numerical value of the similarity parameter with the model nozzle. 

The similarity parameter used in this analysis does not allow an explicit 
determination of tire vehicle base region pressure (P^) and therefore an implicit 
procedure was developed. Figure 3-1 presents a typical base pressure match- 
ing procedure utilizing the base pressure similarity parameter. The illustrated 
curve of possible prototype base pressure ratio as a function of similarity 
parameter can be computed knowing the prototype plume similarity data and 
specifying the flight conditions and corresponding engine chamber pressure 
values. Foi example, Table 3-1 presents prototype plume similarity data 
for a solid rocket booster high performance motor calculated by the use of 
a nozzle flow analysis described in Refs. 2-1 or 2-2. Table 3-2 presents a 
range of flight conditions to be simulated during wind tunnel tests. Choosing 
the flight conditions of P. = 651.4 psia and P = 8.441 psia for M = .796; 

*-'10 oo 

multiplying P^/P by each value of l.O/tP^/P^) of Table 3-1 will yield a proto- 
type possibility curve similar to that of Fig. 3-1. The wind tunnel test is con- 
ducted at the same M = .796 while measurements are obtained at varying 

JO 

model nozzle chamber pressure. Chamber pressure varies directly with 
similar ity parameter . The model nozzle is designed so that the experi- 
mental chamber pressures chosen cause the similarity parameter) , , 

model 

values to fall within the prototype range of interest. The minimum 
value of similarity parameter ) mo( j e j corresponds to that value of P £ that will 
cause flow separation in the model nozzle. This has been experimentally 
determined to occur for P /P^ < 0.60. The maximum value of similarity 
parameter corresponds to either the maximum auxiliary air pressure 
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Similarity Parameter (Varies Directly with P ) 

Fig. 3-1 - Base Pressure Matching Procedure 

1'his is a sample plot and does 
not represent data presented 
in Tables 1-1 and 3 -l 
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Table 3-1 


PROTOTYPE PLUME SIMILARITY DATA FOR A SOLID ROCKET 
BOOSTER HIGH PERFORMANCE MOTOR 


V P b 

.37640+02 
.47970+02 
.61940+02 
.80990+02 
.10732+03 
.14428+03 
. 19694+0 3 
.27 395+03 
.38720+03 
.55923+03 
.82699+0 3 
.12551+04 
. 19609+04 
.31679+04 
.53173+04 
.93368+04 
.17281+05 


M. 

J 

.26599+01 

.27864+01 

.29223+01 

.30761+01 

.32393+01 

.34131+01 

.36005+01 

.38026+01 

.40212+01 

.42585+01 

.45204+0* 

.48153+01 

.51499+01 

.55320+01 

.59676+01 

.64766+0 1 

.70874+01 


6 . 

J 

.10490+02 

.14258+02 

.18025+02 

.21793+02 

.25561+02 

.29329+0 2 

.33096+02 

.36864+02 

.40632+02 

.44395+02 

.48167+02 

.51935+02 

.55703+02 

.59471+02 

.63237+02 

.67006+02 

.70774+02 


.12541+01 

.12685+01 

.12805+01 

.12843+01 

.12887+01 

.12934+01 

.12997+01 

.13062+01 

.13124+01 

.13184+01 

.13248+01 

.13326+01 

.13406+01 

.13485+01 

.13559+01 

.13627+01 

.13687+01 


.17422+02 

.24524+02 

.32211+02 

.40873+02 

.50311+02 

.60603+02 

.71792+02 

.84034+02 

.97486+02 

.11229+03 

.12869+03 

.14695+03 

.16756+03 

.19104+03 

.21793+03 

.24937+03 

.28697+03 


*SP = M. y.) where M = 2.6599. 

.1 J e J e 
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Table 3 -2 

FLIGHT CONDITIONS TO BE SIMULATED 
DURING WIND TUNNEL TESTS 


M 

or 

P./P 

D or 

P 

c 

P 

OP 

SP 

nominal 

.59700+00 

.89000+00 

.72000+03 

.10542+02 

.39132+02 

.79600+00 

.82500+00 

.65140+03 

.84410+01 

.45703+02 

.90000+00 

.77000+00 

.62450+03 

.74250+01 

.50925+02 

.95000+00 

.67500+00 

.61020+03 

.68770+01 

.57365+02 

.10480+01 

.61000+00 

.58000+03 

.57250+01 

.65664+02 

.11000+01 

.64000+00 

.57130+03 

.54770+01 

.64987+02 

. 1 1480+01 

.65000+00 

.57000+03 

.48270+01 

.68890+02 

. 1 2490+0 1 

.69000+00 

.57620+03 

.41050+01 

.73000+02 

.14030+01 

.77000+00 

.58600+03 

.32760+01 

.77946+02 
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supply or maximum mass flow rate capability of the test facility. In the data 
analysis stage of the test program, the model P^/P^ versus similarity param 
ctcr ^model lB P* otte< * on 8arne graph with the prototype possibility curve. 
The intersection of these two curves is the condition which simultaneously 
satisfies both the prototype possibilities and the model actualities and thus 


Similarity Parameter ) model = Similarity Parameter) prototype 


for the given flight conditions. It is this value of P^, P^ that can be expected 
during the actual flight of the prototype vehicle. 


Equation (3.10) presents the general form of the base pressure similarity 
parameter. More specifically, the similarity parameter takes three different 
forms. They are: 


M. 6 


SP 


1 M*~ 5 y. 

<* .1 


J- 


where a = .25 and b = 1.0 


SP 


M. 6. 

2 - 


M* yV 
e 'j 


and 


where a = .25 and b = 0.5 


M. 6. 

SP, = — 


Y i 


where a = 0.0 and b = 1.0 


The actual form of the similarity parameter used in the model nozzle design 
analysis is a function of freestream Mach number and vehicle test configura- 
tion to be considered during the wind tunnel tests. 

Table 3-3 presents a schedule for the recommended use of each form 
of the similarity parameter as a function of freestream Mach number and 
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Table 3-3 

SCHEDULE FOR THE RECOMMENDED USE OF EACH FORM 
OF THE SIMILARI TY PARAMETER AS A FUNCTION 
OF FREESTREAM MACH NUMBER AND VEHICLE 
TEST CONFIGURATION 


M 

aO 

Configuration 

Single Body 
Single Nozzle 

Single Body 
Triple Nozzle 

Triple Body 

0.9 

M. 4. 

C|3 „ J J 

M. 4. 

cn - 4. J , 

M. 6. 
SP = .4 4 

e J 

1 

1 

i.2 

M. 4. 

cp J J 

M. 6. 

cp - -4,4.. 

M. 4. 
SP = . J. ,4 

sp .! - TTUTfroT 

Mo Vj 

S1 1 .,0.25 

M „ y . 
e J 

31 l 6 Ik 

r 

e r j 

1.46 

M. f>. 

SP 2 " .0.2 5 J 0.5 

M e y \ 

M. 6. 

SP 1 ^ 

y j 


3.48 

M 4 

S p ^ — J-J 

3 y i 

M. 6 

SP ^ —1-J. 

3 y j 
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vehicle teat configuration. This table indicates that SP^ is applicable for 
most test cases; SP^ has limited applications; and SP^ is best suited for 
hypersonic Mach numbers. Figure 3-2 defines the variables used in the 
similarity parameters. Figure 3-3 defines the model test configurations 
that may be tested using the recommended similarity parameters. 

When using the computer code to perform the model nozzle design 
analysis, the model nozzles should be designed using the specific similarity 
parameters recommended in Table 3-3. However, it will ususally be found 
that one model nozzle will satisfy all three similarity parameters and there- 
fore the choice of similarity parameter is arbitrary. 

3.3 COMPUTATIONAL STEPS 

Figure 3-1 presents a typical base pressure matching procedure utiliz- 
ing the base pressure similarity parameter. The following paragraphs present 
a detailed discussion of each computational step used to design an envelope 
of candidate model nozzles that may be used in this base pressure matching 
procedure. With the exception of constraints 1 and 2 of Step 1, and the calcu- 
lation of similarity parameters in Step 4, the first four Bteps of the analysis 
are performed by the computer code user and input to the computer code in 
the form of a data deck. (See Section 4-2 for a description of the computer 
code input requirements.) The remaining computational steps have been 
automated by the computer code. 

Step One; 

There are a number of gasdynamic and physical constraints that must 
be considered during the design cf the candidate model nozzles. The gasdy- 
namic constraints are: 

1 . The exit plane Mach number at the nozzle lip (M ) 
must be supersonic. If M„<1 .0, the flow external 
to the nozzle would affect the flow inside the nozzle 
and therefore, the nozzle would not perform as it is 
designed to. And M > 1 .0 allows a simple calcula- 
tion of 6. . e 

3 

3-11 


LOCKHEED HUNTSVILLE RESEARCH & ENGINEERING CENTER 



6 plume boundary initial angle 

plume boundary Maeh number 
M - exit plane Mach number at nozzle lip 
1 . - ratio of specific heats on plume boundary 

M wind tunnel freestream Mach number 

vX 

P wind tunnel freestream statii pressure 

a- 1 

) ratio of specific heats of wine tunnel freestream flow 

jr r 


Fig. l-l - Definition of Variables Used in the Similarity Parameters 
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l. Thu nozzle exit wall angle (6 e ) muni be greater than 
0 deg. Wall angles less than 0 deg can cause strong 
internal shocks resulting in subsonic exit plane Mach 
numbers just discussed. 

3. A gasdynarnic constraint on the ratio P e /P|, is set to 
ensure that model nozzle flow separation does not 
occur. The minimum allowable value of P e /Pb deter- 
mined experimentally is 0.60. 


It is desirable to have one nozzle simulate the flight conditions of a range of 

Mach numbers. This constraint sets a limit on the maximum value of exit 

plane Mach number (M e ) that may be used for the nozzle design. For a 

specified value of P £ and range of flight Mach numbers to be simulated, the 

value of M i is determined at the lowest value of flight Mach number to 
‘max 

be simulated. In most cases, this also corresponds to that flight condition 
where base pressure (P^) is maximum. 

The physical constraints to be set are those imposed by the facility 
where testing is to be performed. 


4. The restrictions of maximum air pressure supply and 
mass flow rate are set. These constraints directly 
affect the maximum obtainable value of similarity 
parameter. This is to say 


Similarity Parameter^ 


max 


oc 


'«), 


oc 


max 




max 


The value of P ) used to determine is 

c /max / max 

either set directly by the test facility or as is most 
often the case, determined as a function of rh) 
as described in Constraint 5 to follow. ^ 


Constraints 3 and 4 are very important considerations in the design of 
the model nozzles. In addition to setting a limit on the maximum value of M 

e 

as already discussed; Constraint 3 also sets a limit on the minimum value of 
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P that may bo used for a given nozzle and set of flight conditions to be 

simulated. Therefore, Constraint 3 sets SP) . in the same manner as 

. /min 

Constraint 4 sets SP) or 

/max 


Similarity Parameter) oc P ) 

/mm c/ min 

Therefore, a given candidate model nozzle may be eliminated from further 
consideration during the analysis simply by comparing the desired similarity 
parameter) prototype to the S]?) min and Sp] max just established by Constraints 
3 and 4 . 


5 . The prototype nozzle exit plane diameter and scale 
of the model nozzle are set. This constraint sets 
the model nozzle exit plane diameter and allows a 
direct calculation of A* for a given model nozzle M 
The value of A“ and the restriction on maximum mass 
flow rate of Constraint 4 yields the value of P \ 

v ymax 

used to calculate SP) described above. P \ 

/max <7 m ax 

is calculated using the isentropic relationship for a 
nozzle flowing air: 



(3.11) 


Step Two: 

Determine a range of flight Mach numbers that is to be simulated during 
testing. Keep in mind that both analytically and physically there is no "one" 
nozzle that may be used to test an infinite range of Mach numbers. Therefore, 
as the range of flight Mach numbers to be simulated increases, the number of 
possible model nozzle solutions decrease and eventually become non-existant. 
This problem, and its solution will be discussed in detail later on in the analysis. 
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Once a range of flight Mach numbers is determined it is necessary to 
choose a schedule of flight Mach numbers to be simulated during testing. 

The schedule chosen will be a function of the test purpose. The actual 
schedule chosen will not affect the final model nozzle design; this has already 
been established when the range of flight Mach numbers to be tested has been 
chosen, Therefore, there is no analytical design limitations introduced by 
having a large number of flight Mach numbers simulated as long as they are 
all within the specified range of flight Mach numbers. 

Step Three: 

Determine the values of freestream static pressure (P ); prototype 

nozzle chamber pressure (P c ); and predicted base pressure ratio (P^/P^ ) 

for each flight Mach number to be tested. The result will be a set of flight 

conditions to be simulated during wind tunnel tests similar to that presented 

in Table 3-2. The values of P and P are easily obtained as a function of 

and/or flight time. Data of this type are presented in motor performance 

prediction and ascent trajectory prediction documentation such as Ref. 3-1 

which is used to construct Table 3-2. The values of predicted base pressure 

ratio used in the analysis are based upon experience gained from previous 

wind tunnel tests. Care should be taken to use the best possible value of 

P^ P for each M . Referring to Step 1, Constraint 3, it is seen that a 

conservative approach may be taken if ratios of P^/P^, slightly higher than 

is actually expected at low are used in the analysis. This approach will 

decrease the maximum allowable value of model nozzle exit plane Mach 

number (M > ) and therefore decrease the size of the family of candidate 

L max 

model nozzles. 


Step Four: 

Determine the prototype plume similarity data for the rocket engine 

that the model nozzle is to simulate. Table 3-1 presents a sample of the 

required data. Figure 3-2 presents a pictorial definition of the variables. 

The tabulated values of M., 6. and y. are the plume boundary properties 

3 1 J 

corresponding to each value of P /P^. These properties are determined 
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utilizing the nozzle flow analysis described in Ref. 2-1 for a gaseous only 
flow, or Ref. 2.2 for a gas. particle flow and an initial plume Prandtl-Meyer 
expansion calculation. The form of the similarity parameter to be used is 
set according to the schedule presented in Table 3. 3. The tabulated value 
of the similarity parameter is calculated by the computer code. 

Step Five: 

Determine the value of similarity parameter required to ensure base 
pressure matching at each scheduled to be simulated. This corresponds 
to ^^* nom j na j * n Table 3-2. First, a curve of possible prototype base pressure 
ratio as a function of similarity parameter is computed as was discussed in 
Section 3.2 and then plotted as in Fig. 3-4. Finally, plotting the predicted 
value of base pressure ratio on the prototype possibility curve yields the 
nominal value of Similarity Parameter (SP ^ required to ensure base 
pressure matching. At this point, it is important to realize that a "predicted" 
value of P^/' P is being used to determine **P nom ^ na j which is in turn used 
to design the model nozzle. Therefore, the design of the model nozzle is only 
as good as the "predicted" value of P^P . This problem and its solution 
will be discussed in more detail later in the analysis. 


Step Five must be repeated for each value of schedule M^ 


lated. Table 3-2 is now complete. 


to be sirnu- 


All data required to begin the model nozzle design compu- 
tations, are now established. The computer code calculation 
is now at what corresponds to Step 7 of Section 2 o£ Point A 
of the computer code flow chart presented in Chart 2-1 . 
These data were read by the computer code in Subroutine 
INPUT. 


Step Six: 

Select: (1) a flight Mach number from the schedule of Mach numbers 
to be simulated as determined in Step 2 and (2) the corresponding flight 
conditions determined in Step 3. 
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The computer code will begin the analysis with the lowest value of 

This will allow an early calculation of M and therefore eliminate the 

1 v. 

max 

needless design calculations for nozzles with an M^ greater than 

e max 

Step Seven: 

Select a value of model nozzle exit plane Mach number slightly greater 
than 1 .0. The value of will be varied parameterically in subsequent itera- 
tions. Selecting a specific value of Mach number gives the model nozzle area 

ratio, A /A . 
e 


Step Eight: 


The next design parameter of concern is the model nozzle chamber 
pressure, P^ . Selecting a specific value of P then permits the calculation 
of plume boundary Mach number according to Eq.(3.6) 



(3.12) 


Note that the pressure on the jet boundary is assumed to be equal to the base 

pressure. P. is determined as a function of the predicted value of P. / P 
r b r o oo 

in Step 3 and P^) , . in the wind tunnel. Please note that tire P_) 

model 00 prototype 

is not used in this calculation. The value of P £ chosen in this step is held 

constant throughout the analysis and is usually set equal to 50% of the maximum 

available air pressure supply. This value of P is chosen to ensure that each 

model nozzle in the family of candidate nozzles being designed will be capable 

of producing a power sweep curve with a range of similarity parameter values 

bracketing SP . (See Fig. 3-1). The value of P is read in Subroutine 

nominal n c 

INPUT. 
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Step Nine: 

The value of 6. can be found which ensures similarity using Eq.(3.9) 


A. 

J 


^ model 


- SP 


nominal 


) 


prototype 


^ 1 model 


(3.13) 


where: ^ = 1 .4 for a model nozzle flowing air; a and b are functions of flight 

Mach number and vehicle test configuration and are determined in accordance 
with the schedule presented in Table 3-3. This value of 6^ will ensure base 
pressure matching at the specified set of flight conditions, the chosen value 
of and the constrained value of P only . 


Step Ten: 

Finally, the model nozzle exit plane wall angle, 0 , required to achieve 

6 ) for the specified conditions is calculated according to Eq.(3.8) 

model 


0 

e 


• model 


where u and 
e 


are calculated according to Eq. (3.5) 


and 




- tan 



1 ) 


1/2 


"j 



tan 



1) 


1/2 


- tan" 1 (M? - l) 1 ' 2 
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At this point, "a" candidate nozzle with "a" defined value of 0 and M and 

e c 

"a" constrained value of P £ has been designed which will ensure base pres- 
sure matching for " one" specified set of flight conditions. Should one flight 
condition change, the result would be a change in ®P nom j na j ant * therefore 
base pressure matching would not occur. To regain the condition of base 
pressure matching requires only a change in P c proportional to the change 
in SP . . This assumes that the "new" SP . does satisfy the 

condition . 


SP . . < 

minimum — 


new 


SP . < SP 

nominal — maximum 


for the given candidate nozzle. 

Step Eleven: 

Incrementing the value of and repeating Steps 9 and 10 will 
define a family of candidate nozzles which will: (1 ) ensure base pressure 
matching for "one" specified set of flight conditions; and ( 2 ) meet all of the 
constraints specified in Step 1 . Figure 3-5 presents a sample plot of a 
family of candidate model nozzles that may be used to ensure base pressure 
matching for the flight conditions to be simulated for M^ - .597 present in 
Table 3-2. Note that nozzle flow separation occurs for greater than 3.83. 

Step Twelve: 

Repeat Steps 6 through 1 1 for each flight Mach number in the schedule 
of Mach numbers to be simulated. The result will be one family of candidate 
model nozzles that may be used to ensure base pressure matching for each 
of the specified set of flight conditions to be simulated and can be plotted 
similar to that in Fig. 3-5. 

Step Thirteen; 

Each family of candidate model nozzles derived for a particular flight 
Mach number is now checked to see which nozzles in that family are capable 
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maximum 


5.0 


3.6 3.8 4.0 


y be Used to Ensure Base 
= .597, Pj/P^ = .89, 

>2 osia and SP . = 39. 1 32 



LMSC-HREC TR D784111 


of simulating the flight conditions that exist for each of the other flight Mach 
numbers. This is achieved as follows: 


Select a family of candidate nozzles derived in Step 12. For each 

candidate nozzle in that family, calculate the values of SP . . and 

7 minimum 

SP that are obtained when tested at each of the other fligh Mach 

maximum 

numbers. Then, compare these values to the required SP , that corre 

r ^ nominal 

sponds to each flight Mach number to ensure that the constraint 


SP . + 

minimum 


5.0 < SP 


or . , 

— nominal 


SP 




5.0 


(3.14) 


is met (see Fig. 3-1 for a clear understanding of the importance of this 
const raint) . 


The result is a group of nozzle families that form an "envelope of model 
nozzles" that may be used to ensure base pressure matching for each set of 
flight conditions to be simulated. There is one envelope of nozzles for each 
set of flight conditions to be simulated. 


The tolerance of + 5.0 on SP and -5.0 on SP in Eq.(3.l4) 

minimum maximum n 

is applied to introduce a little more conservatism in the final nozzle design. 
These tolerances will: (1) allow for changes in SP , that may result if 

wind tunnel test conditions vary slightly, and (2) decrease the accuracy re- 
quired on the initial predicted values of P^ P used in Step 3. 


Refer to Fig. 3-5 and note that each candidate nozzle represented along 
the curve will ensure base pressure matching for SP , = 39.1 32. How- 

ever, when the above tolerances are applied, the final accepted nozzles in this 
family for this set of flight conditions become those that are bracketed by 

UZ 111 • 
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The values of SP . . and SP # _ . are calculated using the 

minimum maximum 

general form of the Similarity Parameter 


SP) 


M. 6. 


general ~ M a b 

e j 


For a given value of and a nozzle flowing air y^ = 1.4; the minimum 
and maximum values of S p )g enera i are determined l- imply by determining the 
minimum and maximum values of M. 6. respectively. 


j J 


Determination of SP) 


minimum 


a. M , 0 , P P and v are defined for each candidate nozzle, 
c e c e e 

P 


c \ 

>. M.\ and 6.\ occur for 1 

^'min -Vmin b / 

' min 


c . 


For a given P /P^; (P / P (} ) can be determined using 
constraint 3 of Step 1 


min 


P 


run 


> \ P 

S) = _£ 0 

) I T-. V. 

b/ . 


mm 


For a given set of flight conditions, P^ is given and P \ 

/ m i 

can be calculated from the above equation 


min 


c) . “ ^ * 

'min « 


0.60 t P, 
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d. 


M.\ and 6.) arc then determined using Eqs.(3.6) 

^'min ^*min 

and (3.7), respectively 




- I 


1/2 



0 

e 



min 


v 


e 


where 


p.\ is th ter mined using Eq.(3.5) 
•*'min 



finally, 


e. SP . 

minimum 


m j) 


min 




. ,a ET 
M y . 
e J 


min 


* 


Determination 



maximum 


P . 

a. M , G • "pr- arid p is defined for each candidate nozzle, 
e e P e 


b. 


M i) 


max 


and 


6 . 

J 


max 


occur 


for 



max 
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c. For a given flight condition is defined. Therefore 
P £ /P^ is maximum when P c is maximum. The value 

of P ] is the lesser of the test facility Pi or 
c * max i c / max 

as ts most often the case, the value of P determined 

. \ c ' max 

as a function of ml for the facility 

/max 7 




nil 1 

Jt~ 

/ max 

1 o 


Therefore , 


/ max 

P L 


d. M 


M ) , i< \ and 6 \ are determined using Eqs, (3.6), 

'ma x *’ 'ma x ^ ' max 

( >.5), and (3.7), respectively. 




V: * 1 


V: 


i \ = 0 + i >. ) - v 

i / e j / c 

•"max "max 




■,) ■{&) 




- tan 


- 1 M 2 ) 

j L 
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Finally, 


e . 


SP 


maximum 


m j) * *j) 

J f max J > ma x 

M a y h ~~ 
c r ) 


Step Fourteen 

The results of Step 13 are plotted. Figure 3-6 presents the results of 
Step 13. For example. Fig. 3-6a presents an envelope of candidate model 
nozzluB that may be u^ed to simulate the flight conditions that exist for M^ 

- 0.597. Each candidate nozzle in this envelope will: (1) ensure base pres- 
sure matching for the = 0.597 flight conditions, and (2) meet all of the 
constraints of Step 1 . 


Steps 15 through 18 compare each of the envelopes of nozzles plotted in Fig. 
3-6 to determine a "final" envelope of model nozzles that may be used for 
all flight Mach numbers to be simulated. 


Step Fifteen 


The upper limit on is a result of the separation constraint of Step 1 

and visually occurs for the lowest value of flight Mach number to be simulated. 

Compare each envelope of model nozzles and choose the smallest value of M^ 

as the upper limit on M^ of the "final" model nozzle envelope. This would 

correspond to M - 3.53 in Fig.3-6a for M - .597. 
r e oo 


-max 


Step Sixteen 

The lower limit on M^ is a result of the mass flow constraint of Step 1 . 
This constraint has its greatest influence at the highest value of flight Mach 
number to be simulated. This can be understood by examination of Table 3-2. 
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A« flight Mat h number to be simulated increases, the value of SP nQ 

required to ensure base presHure matching also increases. P £ is directly 

proportional to ®P nom [ na i an< * therefore also must increase with increasing 

flight Mach number to be simulated. For a given P £ and constrained nozzle 

exit diameter the mass flow rate increases as M decreases. 

c 


Therefore, examine the nozzle envelope for the highest flight Mach 

number to be simulated. Compare each of the nozzle families and choose 

the largest value of M as the lower limit on M of the "final" model 

e min 

nozzle envelope. This would correspond to - 3.49 in Fig. 3.6i for = 

1.403. Once the upper and lower limits of model nozzle exit Mach number 

are determined the value of the limits are compared. If the lower limit on 

M is greater than the upper limit on M p no solution exists for the range of 

flight Mach numbers to be simulated (see step two). To solve this problem 

merely reduce the upper value on the range ?f flight Mach numbers to be 

simulated during testing established in Step 2. This will reduce the value of 

M just determined to some value less than M of Step 15. 

e min e max 


In some instances, as is the case of the example being presented, the 

lower and upper limits of M are close to the same value (3.49 and 3.53). 

c 

Further reduction in the upper value on the range of flight Mach number to 
be simulated during testing will increase the difference of M o and M^ 
of the "final" model nozzle envelope. 


'min 


max 


Step Seventeen 

Determine the maximum and minimum values of 0 for all candidate 

e 

nozzles whose M g correspond to the upper limit of M g determined in Step 15. 
Step Eighteen 

Determine the maximum and minimum values of 0 for all candidate 

e 

nozzles whose M g correspond to the lower limit of M fi determined in Step 16. 
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NOZZL. E FAMILY 


MBOL 

I DENT I F I c AT ! 

ON 



□ 

DERIVED 

FOR 

FLIGHT 

MACH 

NO. 

= 

0.597 

A 

DERI VL'J 

FOR 

flight 

MACH 

NO. 

3 

0.706 

0 

DERIVED 

FDR 

FLIGHT 

MACH 

NO. 

E 

0.90 

0 

DERIVED 

FOR 

FLIGHT 

MALH 

NO. 

= 

0.95 

A 

DERIVED 

FOR 

flight 

MACH 

NO. 

3 

1 . 048 

L 

DERI VED 

FUR 

Flight 

MACH 

NO. 

s 

1.10 

# 

0! RIVED 

FUR 

flight 

MACH 

NO. 

a 

1.148 

V 

D! RIVED 

FOR 

flight 

MACH 

NO. 

3 

1.249 

0 

DERIVED 

FOR 

flight 

MACH 

NO. 

= 

1.403 



3 .>* 3 . 

EXIT MACH NUMBER 


p ig. 3-6a - The Symbols Plotted Above Represent Model Nozzles that May 
be Used to Simulate Only the Flight Conditions that Exist for 
a Freestream Mach Number = .597 
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NOZZLE FAMILY 
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Fig. 3- 6b - The Symbols Plotted Abo/e Represent Model Nozzles that May 
be Used to Simulate Only the Flight Conditions that Exist for 
a Freestream Mach Number =•. .796 


3-30 


KVKHitn HllNISVILtl HI '.t AHl M & t NGINt l KING CfcNTl R 







Fig. 3-6c - The Symbols Plotted Above Represent Model Nozzles that May 
be Used to Simulate Only the Flight Conditions that Exist for 
a p i eestream Mach Number = .900 
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Fig. 3- 6c - The Symbols Plotted Above Represent Model Nozzles that May 
be Used to Simulate Only the Flight Conditions that Exist for 
a Freestream Mach Number = .950 
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Fig. 3-6f - The Symbols Plotted Above Represent Model Nozzles that May 
be Used to Simulate Only the Flight Conditions that Exist for 
a Freestream Mach Number = 1.100 
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Fig. 3-6g - The Symbols Plotted Above Represent Model Nozzles that May 
be Used to Simulate Only the Flight Conditions that Exist for 
a Freestream Mach Number = 1.148 
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Fig. 3-6h - The Symbols Plotted Above Represent Model Nozzles that May 
be Used to Simulate Only the Flight Conditions that Exist for 
a Freestrcam Mach Number = 1.249 
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Fig. 3-6i - The Symbols Plotted Above Represent Model Nozzles that May 
be Used to Simulate Only the Flight Conditions that Exist for 
a Freestream Mach Number = 1.403 
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Step Nineteen 

The results of Steps 15 through 18 can be plotted to yield the final 
envelope of model nozzles that may be used in the test program. Each model 
nozzle in this envelope wills (1) ensure base pressure matching for "all" 
specified sets of flight conditions to be simulated, and (2) meet all of the 
constraints specified in Step 1. Figure 3-7 presents a plot of the final envelope 
of nozzles that result, when Steps 15 through 18 are applied to Figs. 3-6a through 
3-6i. 


Step Twenty 

To ensure that base pressure matching docs occur in the data analysis 

stage of the test program, it is necessary to know the power sweep operating 

characteristics of the model nozzle prior to actual testing. For a given 

model nozzle design, the power sweep operating characteristics of importance 

are: (!) the minimum obtainable value of similarity parameter (SP . . ); 

(2) the maximum obtainable value of similarity parameter (SP ) for 

each of the flight conditions to be simulated; (3) t lie nominal value of similarity 

parameter (SP , ,) required to ensure base pressure matching and (4) the 

nominal ° 

values of chamber pressure corresponding to each value of similarity param- 
eter just noted. 


The calculations of SP 


maximum 


, p c ) 

’ 'max 


, SP . . and P 
' minimum c 


/mi 


have 


min 


already been discussed in Step 13. The value of SP nom ^ na j required to ensure 

base pressure matching has already been defined in Step 5 for each scheduled 

M to be simulated. This leaves the value of chamber pressure, P J , 

^nom 

corresponding to ®P nom ^ na j as the °nly value undefined at this time. 


Refer to Fig. 3-1 when performing this step of the analysis. 
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The value of P \ is calculated following the interation process 
mom 

is a function of flight conditions to be simulated 

nom 

as well as model nozzle geometry. 

P c * 

a. M , 0 , , D , D and i; are defined for each candidate 

C 6 tr © © 

© 

nozzle. 

b. Set a value for the initial guess on M. 

J 

c. Calculate 6. using Eq.(3.10) 

J 

M* y b 

6 j “ M. J * ® P nominal 
J J 


d. Calculate 0^ using Eq. (3.8) 


0 = 6. + v - j/. 

c J e j 


where v j is determined using Eq. (3.5) 



e. Compare the value of 0 C just calculated with the actual value 
of 0 e for the model nozzle being considered. If the values of 
0 e compare within an acceptable tolerance proceed to f; if 
not, adjust the value of M: and repeat steps c, d and e until 
an acceptable tolerance is achieved. 

P 

f. At this point M: is defined and the value of p- is determined 

using Eq. (3.3) J b 
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g. Finally, P.i is calculated using tlu* following relationship 
'nom 



nom 



❖ p \ 

80 ' model 


Table 3-4 presents power sweep operating characteristics of five 
representative model nozzles that exist in the final envelope of Fig. 3-7. 
The flight conditions fro be simulated during wind tunnel tests are those of 
Table 3-2. 


Examination of Table 3-4 indicates that each of the representative 
model nozzles considered are capable of base pressure matching well within 
the constraints of Step 1. This is indicated by the fact that the constraint 


SP 


minimum 


+ 5.0 < SP 


< SP 

nominal — maximum 


- 5.0 


is satisfied for all model nozzles at all scheduled M to be simulated. 

00 

The analytical procedure used to design model nozzles which meet 
MSFC base pressure similarity parameter criteria is now complete. The 
following section presents an input guide for an automated computer code 
used to perform the described analytical procedure. It is strongly suggested 
that the computer code user make every possible use of the analytical portion 
of this document when setting up his input deck. Doing so will give the user 
a greater understanding of model nozzle design problem and the actual use 
of the computer code. 
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4. A COMPUTER CODE USED TO DESIGN MODEL NOZZLES 
WHICH MEET MSFC BASE PRESSURE 
SIMILARITY PARAMETER CRITERIA 

This section describes a computer code used to design model nozzles 
which may be used to simulate prototype engine conditions that exist over a 
specified range of ascent trajectory conditions. Each of the model nozzle 
designs will be capable of achieving the desired similarity parameter values 
required to ensure base pressure matching within the restrictions of maximum 
air pressure supply and mass flow rate imposed by the facility where wind 
tunnel testing will be performed. 

An example of the computer code input requirements and resultant out- 
puts are presented for the wind tunnel test LA-604 planned for the NASA-Ames 
Research Center 11x11 wind tunnel (Ref. 1-1). A "final envelope" of model 
nozzle designs, and corresponding power sweep operating characteristics for 
each of the desired flight conditions to be simulated in the wind tunnel tests 
is presented for representative model nozzles in the "final envelope." 

Subroutines that make up the computer code are listed in Table 4-1. 

The subroutines which call and are called by the particular routine as well 
as a brief statement regarding the function of the routine are also included 
in the table. Chart 4-1 presents a flow chart of the computer code. A listing 
of the computer code is presented in the Appendix. 

Sections 2 and 3 of this report present a detailed development of analytical 
procedures used in the model nozzle design analysis and automated in this com- 
puter code. It is strongly suggested that the computer code user make every 
possible use of the analytical portion of this report when performing a model 
nozzle design analysis. By doing so, the user will gain a greater understanding 
of: (1) the model nozzle design problem being considered, and (2) the actual 
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Provides overall program 
control 


Data input control and 
definition of model 
nozzle design problem 


Controls derivation of each 
"family" and then each 
"envelope" of candidate 
model nozzles 


Controls derivation and 
plotting of the "final 
envelope" of candidate 
model nozzles and deter- 
mination of power sweep 
operating characteristics . 


Controls plotting 
of optional plots 


r r — 

I | MORE 


I 


DESIGN 


J 



Chart 4-1 - Flow Chart of a Computer Code Used to Design Model 
Nozzles Which Meet MSFC Base Pressure Similarity 
Parameter Criteria 
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use of the computer code. This section of the report will make continuous 
reference to the applicable paragraphs of Sections 2 and 3; this will allow the 
user to go directly to the paragraphs which will help answer possible questions. 

4.1 COMPUTER CODE CAPABILITIES 

The computer code described in this section automates the analytical 
procedures described in Sections 2 and 3 and is used to design an "envelope’' 
of model nozzles (the envelope is defined by model nozzle exit plane Mach 

number, M , and model nozzle exit wall angle, 6 ) that may be used in a 

c e 

specified wind tunnel test program. Each model nozzle in this envelope will: 

(1) be capable of achieving the desired similarity parameter values required 
to ensure base pressure matching for all specified sets of flight conditions to 
be simulated, and (2) meet all the gasdynamic and physical constraints that 
must be considered during the design of a model nozzle. (See Section 3.3, 

Step One, for a discussion of these constraints.) 

The analytical procedure discussed in Sections 2 and 3 is not limited to 
Space Shuttle vehicle configuration applications. Consequently, this computer 
code may be used to design model nozzles to simulate prototype engine condi- 
tions of any single body-single nozzle, single body-triple nozzle or triple-body 
engine system. (See Fig. 3-3 for a definition of these engine systems.) 

To ensure that base pressure matching occurs in the data analysis stage 
of the test program, it is necessary to know the power sweep operating char- 
acteristics of the model nozzle prior to testing. The computer code will calcu- 
late the power sweep operating characteristics of representative model nozzles 
that exist in the "final envelope" of model nozzle designs for each of the desired 
flight conditions to be simulated during the test pregram. Section 3.3, Step 
Twenty, lists the power sweep operating characteristics of importance, presents 
their method of calculation, and shows how they may be used to choose a final 
model nozzle design that will ensure base pressure matching at all flight condi- 
tions to be simulated. 
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The final model nozzle design chosen is defined by a model nozzle exit 

plane Mach number, M , and a model nozzle exit wall angle, 6.. The wall 

c c 

contour of the actual model nozzle used in the wind tunnel test is arbitrary 
and left up to the nozzle designer. Within the scope of the current technology 
it is felt that a conical nozzle would suffice for all cases. 
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4.2 USER'S INPUT GUIDE FOR THE MODEL NOZZLE DESIGN COMPUTER 
CODE 

This section outlines in detail the procedures for using the model nozzle 
design computer code. Each card and its use are explained in Section 4.2.1 . 
The control card set-up for the Univac 1108 Exec 8 is presented in Section 
4.2.2. 

4.2.1 Computer Code Input Information 

The input data are organized into sections determined by their use. 

An example of the computer code input requirements used to design an 
"envelope" of model nozzles for possible use in wind tunnel test IA-604 
planned for NASA- Ames Research Center 11x11 wind tunnel (Ref. 1-1) is 
presented in Table 4-2. A description of these cards is given below. 


Computer Code Input Instructions 


Card 1 Run Control Card 

Column Parameter Value 
5 NMACH 


10 NSP 


1 


Format 515 (Right Adjusted) 

Description 

Refer to Section 3.3, Steps 2 and 16. Antici- 
pated number of scheduled flight Mach numbers 
to be simulated during testing. If the range 
of flight Mach numbers to be simulated is 
such that no solution exists the upper value 
on the range of flight Mach numbers will be 
reduced automatically until a nozzle design 
solution does exist . The value of NMACH 
will be adjusted accordingly internal to the 
computer code (20 max.). 

Refer to Section 3.2. Controls the form of 
the similarity parameter used in each nozzle 
design analysis. The similarity parameter 
to be used in the analysis is specified accord- 
ing to the schedule recommended in Table 3-3. 

M. 6. 

SP = — sU is used in the analysis 

M 0 * 25 y. 
e J 
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Cak'd 1 (Cont'd) 

Column Parameter Value 


10 NSP 2 


3 

1 5 NPRINT 0 


1 

20 NPEOT 0 


1 


2 5 NNOZF 1,2,3, 

4 or 5 


Description 

M. 6. 

SP - -- - is vised in the analysis 
M 0.« 0.5 

« J 

M. 6. 

SP = — is used in the analysis 


Refer to Section 3.3, Step 20. There will 
be no intermediate printout. Only the power 
sweep operating characteristics of repre- 
sentative model nozzles for each of the flight 
conditions to be simulated will be printed. 

Print intermediate data showing the results 
of the calculations in Section 3.3, Steps 12 
and 14 (Fig. 3-7). 

Refer to Section 3.3, Steps 15-19. There 
will be no intermediate plots. Only the 
"final envelope" of model nozzles that may 
be used in the test program will be plotted. 
Each model nozzle in this envelope will; 

(1) ensure base pressure matching for all 
specified sets of flight conditions to be 
simulated, and (2) meet all the constraints 
specified in Section 3.3, Step 1 . 

Plot intermediate data showing the results 
of the calculations in Section 3.3, Step 1‘* 

(Fig. 3-7). 

Number of representative model nozzles in 
Fig. 3-8 that power sweep operating char- 
acteristics will be determined for (max. of 
5). 

Example; If NNOZF - 4, power sweep opera- 
ting characteristics will be determined for 
Nozzles 1, 2, 3 and 4. 
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Cards 2 Axes Identification Labels for the "Final 
Envelope" of Model Nozzles 

The following set of cards contains the labels placed on the horizontal 
and vertical axes when the "final envelope" of model nozzles that may be used 
in tht test program is plotted. Please note that model nozzle exit Mach number 
and model nozzle lip angle are always plotted on the horizontal and vertical 
axes, respectively. Cards 2 merely allow the program user to clearly identify 
the axes labels with the model nozzle design problem being considered. (Refer 


to Fig. 4- 

1). 



Card 

Column Parameter 

Description 

Format 

2.1 

1-48 XFINAL 

Horizontal axis label 

8A6 

2.2 

1-48 Y FINAL 

Vertical axis label 

8A6 

Cards 3 

Vertical axis and nozzle design problem identification for optional 
plots (The following cards are not required if NPLOT = 0). 


The following set of cards contain the labels placed on the vertical axis 
on each of the optional plots and problem identification information. These 
labels do not change from plot to plot. As in Cards 2, the nozzle lip angle is 
always plotted on the vertical axis and Card 3.2 merely shows the program 
user to more clearly identify the axes labels with the model nozzle design 
problem being considered (Refer to Fig. 4-2 and 4-2a). 


Card 

Column 

Parameter 

Description 

Format 

3.1 

1-48 

TITLE 

Nozzle design problem 
identification label 

8A6 

3.2 

1-36 

XLINE 

Vertical axis label 

6A6 


Cards 4, 5 and 6 contain the prototype plume similarity data for the 
rocket engine that the model nozzle is to simulate. These data are calculated 
in Section 3.3, Step Four. Table 3-1 presents a sample of the required data 
(the column labeled SP is not required input). Figure 3-2 presents a pictorial 
definition of the parameters input on these cards. The values of M., 6. and y. 
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are the plume boundary properties corresponding to each value of P^yP^. 
These properties are determined by means of a numerical integration 
through a Prandtl- Meyer expansion. 


Card 4 - Format El 0.5 
Column Parameter 

l - 1 0 PXME 

Card 5 - Format 15 (Right 
Column Parameter 

1-5 NPMRAY 


Description 

prototype exit plane Mach number at nozzle 

lip. 

Adjusted) 

Description 

Number of plume boundaries for which plume 
similarity data are given (30 max.) 


Cards 6 Prototype Plume Similarity - Format 4E10.5 

Data (each card contains the required similarity 
data along one plume boundary and is repeated 
NPMRAY times). Card 6 is input in order of 
increasing value of PPCOPB. 


Column Parameter Description 

1-10 PPCOPB P^P^ " P rot:0 tyP e pressure ratio where P £ is 

nozzle chamber pressure and P^ is the pressure 
along the plume boundary, P^P^ may be con- 
sidered the independent variable, with the remain- 
ing parameters on this card determined as a func- 
tion of P /P,. 

c b 


1 l -iO 
<il -30 
31 -40 


PXMJ 

PDELJ 

PGAMAJ 


- Prototype plume boundary Mach number 
corresponding to P /P^. 

6^ - Prototype plume boundary initial expansion 
angle corresponding to P c /P^. 

Tj - Prototype ratio of specific heats on the 
plume boundary corresponding to P /P^. 
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Cards 7 and 8 contain both the program constants and the gasdynamic 
and physical constraints that must be considered during the design of the 
model nozzle. The constraints that are input are described as Constraints 
3,4 and 5 in Section 3.3, Step Two. Constraints 1 (M > 1 .0) and 2 (6 > 0) of 
Step 2 arc set internal to the computer code and are not input requirements. 


Card 7 - Format 7E10. 5 
Column Pa ramete r 

1.10 PCMAX 


11-20 PC 


Description 

P ) - Refer to Section 3.3, c tep 1, Constraints 

c 'max 

4 and 5. Maxi' i o model nozzle chamber pres- 
sure (psia). This corresponds to the maximum 
air pressure aivf'iy capability imposed by the 
facility where testing is to be performed. This 
value serves as an upper limit only. The actual 
maximum value of P is calculated as a function 

of A*, T and m\ _ . 

’ o /max 

P c - Refer to Section 3.3, Step 8. This is a 
specified value of model nozzle chamber pressure 
(psia). This value is held constant throughout 
the nozzle design analysis and is usually set at 
0.50 * PCMAX. 


21-30 


31-40 


PEPBMN 


PATM 


P /P, ) - Refer to Section 3,3, Step 1, 

c D /min 

Constraint 3. Minimum allowable value of 
P e /P^ set to ensure that nozzle flow separation 

does not occur. The value of P /P, ) deter- 

e/ b/ min 

mined experimentally is 0.60. 

Wind tunnel total pressure (psia). 


Note : The nozzle design analysis is performed assuming that the model 

nozzle will be flowing a constant }' simulant gas during wind tunnel 
tests. Therefore, the values of GAMAJ (Column 41-50) and GAMAE 
(Columns 51-60) must be equal and are held constant throughout the 
nozzle design analysis. 


41-50 


GAMAJ 


Tj - Refer to Fig. 3-2. Model nozzle ratio of 

specific heats on the plume boundary. For a 
nozzle flowing air y. - 1 .40. 
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Column Parameter 
51-60 GAMAE 

61-70 GAMA1 


Description 

y - Refer to Fig. 3-2. Model nozzle ratio of 

specific heats at the nozzle exit plane. For a 
nozzle flowing air = 1 .40. 

Too - Refer to Fig. 3-2. Ratio of specific heats 
of wind tunnel freestream flow. For air y ^ = 1 .40. 


Format 4E1 0.5 
Parameter 


Card 8 - 
Column 
l -10 

11-20 

21-30 

31 -40 


SCALE 

PEX1TD 

WMAX 

TO 


Description 

Refer to Section 3.3, Step One, Constraint 5. Ratio 
of model nozzle exit diameter to prototype nozzle 
exit diameter. This constraint sets the model 
nozzle exit diameter and allows a direct calculation 
of A* for a given model nozzle M e . 

Prototype nozzle exit diameter used to define 
SCALE (in.) 

m) - Refer to Section 3.3, Step One, Constraint 
'max 

4. Maximum model nozzle mass flow rate. This 
corresponds to the maximum air mass flow rate 
capability imposed by the facility where testing 
is to be performed (lbm/sec). 

T - Refer to Section 3.3, Step 1, Constraint 5. 
o 

Model nozzle chamber temperature used in the 
P ) calculation (R). 
c max 


Card 9 Flight Simulation Data and Labels for Optional 
Plots (the following cards are repeated NMACH 
times). Cards 9 are input in order of increasing 
value of XMI. 


The following set of cards contain the values of freestream static pressure 

(P ), prototype nozzle chamber pressure (P ), and predicted base pressure 
oo c 

ratio (P. /P ) determined in Step 3 for each flight Mach number to be simulated 

o oO 

during testing as established in Step 2. The result will be a set of flight conditions 
to be simulated during wind tunnel tests similar to that presented in Table 3-2. 

If NPLOT equals 1, additional labels for the optional plots are also input at this 
time. 
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Car d 9. 1 
Column 
1-10 

11-20 

21-30 

31-40 

Card 9 .2 

Column 

1-60 

Card 9.3 

Column 

1-36 


Flight Simulation Data 
Parameter 


XMI 


PBOPI 


- Format 4EI 0.5 
Description 

M - Refer to Section 3-3, Step 2, and Fig. 3-2. 
00 

Scheduled flight Mach number to be simulated 
during the wind tunnel test. 

P^/P^ - Refer to Section 3-3, Step 3. Predicted 

base pressure ratio to occur at the specified 
value of XMI in the ascent trajectory. 


P c - Refer to Section 3-3, Step 3. Prototype 

nozzle chamber pressure which occurs at the 
specified value of XMI in the ascent trajectory 
(psia). 

P - Refer to Section 3-3, Step 3, and Fig. 3-2. 

uO 

Prototype freestream static pressure which 
occurs at the specified value of XMI is the ascent 
trajectory (psia). 

- Nozzle Family Identification Label - Format 10A6 
(not required if NPLOT equals 0) 

Parameter 

NOZD Refer to Fig. 4-2 and Section 3-3, Step 11. This 

label identifies the family of candidate nozzles 
which have been designed to ensure base pressure 
matching for the flight conditions that exist at the 
specified value of XMI on Card 9.1 . 


PPC 


PPI 


- Horizontal Axis Label - Format 6A6 
(not required if NPLOT equals 0) 

Parameter 

HEDRP Refer to Fig.4-2a. Horizontal axis label for 

the specified value of XMI on Card 9.1 . Please 
note that model nozzle exit Mach number is a 
always plotted on the horizontal axis. Card 9.3 
merely allows the program user to more clearly 
identify the axis label with the model design 
problem being considered. 
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Table 4-2 


AN EXAMPLE CASE SHOWING THE REQUIRED INPUT FORMAT USED 
TO DESIGN AN "ENVELOPE" OF MODEL NOZZLES FOR POSSIBLE 
USE IN WIND TUNNEL TEST IA-604 PLANNED FOR NASA-AMES 
RESEARCH CENTER 11*11 WIND TUNNEL (REF.l-t) 


Card 1 
Card 2.1 
Card 2.2 
Card 3.1 
Card 3.2 
Card -1 
Card 5 


Card 6 
Repeated 
NPMRAY 
Time* 


Card 7 
Card 8 
Card 9. 1 
Card 9.2 
Card 9.3 


Repeated 

NMACH 

Times 



9 1 

1 1 

5 


.02 SCALE 

HPM nozzle 

EXIT MACH NUMBEI 

.02 SCALL 

HPM nozzle 

LIP ANGLC 

.02 SCALE 

HPM NOZZLE 

0ES1GN 

ANAL VS1S 

NOZZLE L1R 

' ancle 



2.6599 




17 




37.64 

2.6599 

l0.9«0 

1 .2541 

47.97 

2.7864 

14.258 

1.2685 

61.94 

2 .9 22 3 

18.025 

1.2805 

80. 99 

3 . 0 76 1 

21.793 

1 .2843 

107.32 

3.2 39 3 

25.561 

1.2807 

144 .28 

3.4 131 

29.329 

1 . 29 34 

196 .94 

3.6005 

33.096 

1.2997 

273.55 

3.8 026 

36.864 

1 . 3062 

387 .2o 

4.C 212 

40.632 

1 . 3124 

559.23 

4.2 585 

44.395 

1 . 3184 

8 26 .99 

4.5 204 

48.167 

1 . 32 48 

1255.1 

4.8 153 

51.935 

1 . 3326 

196C.9 

5.1 499 

55.703 

1 . 3406 

3167.9 

5.5 320 

59.471 

1 . 3485 

531 7. ? 

5.9676 

63.237 

1 . 3559 

9336.8 

6.4 766 

67.006 

1 . 3627 

17281.0 

7.0e74 

70. 774 

1 . 3687 

15UT.C 

750 .0 

0.60 

14.70 

.02 

149 .64 

27. ci 

50 C. a 

.597 

.89 

720.0 

10.542 


TEST IA-604 


no? ? i c fahily oERivro for flight mach 

N027LE E * I T MACH NUHBEr MINF=0.597 
•796 >62! 651.* 6.941 

NOZZLE FAMILY OEBIVFD FOR FLIGHT MACh 
N0Z7LE E * I T MACH NUMBER MINF:0.796 
.90 .77 62 *. 5 7.42b 

N027LE FAMILY OERIVCO FOR FLIGHT MACH 
NO? ?LE EXIT MACH NUMBER MINF:0.90 
•95 .675 610.2 6.877 

N 0 ??LC FAMILY DERIvrO FOR FLIGHT MACH 
N02 ?Lf EXIT m*CH NUMBER MINF=0.95 


1 .40 


NO. - 0.597 


NO. = 0.796 


NO , 


0.90 


NO. : 0.95 


1.04b .61 580.0 

N022CE FAMILY DERIVED FOR 
N022LE EXIT MACH NUMBER 
1.10 .64 S 7 1 . 3 

N022LE FAMILY OERIVED FOR 
NOZ ?LC EXIT MACH NUMBER 
1.148 .65 570.0 

N0??LE FAMILY OERIvEO FOR 
NOZZLE EXIT * *CH NUMBER 
1.249 .69 576.2 

NOZZLE FAMILY DERIVED FOR 
NOZZLE EXIT MACH NUMBER 
1.4C3 .77 586.0 

NOZZLE FAMILY OERIVEO FOR 
NOZZLE EXIT MACH NUMBER 


5. 725 

FLIGHT MACh NO. 
M INF : 1 .04 8 
5.477 

FLIGHT MACH NO. 
MINF = 1 .10 
4.827 

FLIGHT MACH NO. 
M INF £1.148 
4.105 

FLIGHT MACh NO. 
MlNF : 1 .249 
3.276 

FLIGHT MACH NO. 
M INF : 1 ,4 LJ 


1.048 
1.10 
1.148 
1 .249 
1 .403 


l.*»0 


1.40 
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4.2.2 Control Card Set-Up for Univac 1108 Exec 8 

A typical run stream set-up for the Univac 1108 Exec 8 computer is 
presented in this section. It is presented to acquaint the program user with 
the magnetic tape assignments and required communications with the com- 
puter. The data deck has been described in Section 4.2.1 and an example 
case presented in Table 4-2. 

The program tape referenced on the following page contains a data file 
designated ae PPM. The data in this file are presented in T • ble 4-2 for the 
wind tunnel test IA-604 planned for NASA-Amcs Research Center 11 x 11 ft 
wind tunnel (Ref. 1-1). Therefore, if the program user wishes to run a sample 
case, all that is necessary is to replace the DATA DECK cards with the one 
card: 


7 ADD, P HPM 


If the program user desires to add his own new data file to the progi .n 
tape the following cards must be inserted as indicated by "optional insert" in 
the runstream. 

( 7 END 

"Optional Insert" ( NEW DATA DECK (according to Section 4.2.1) 

( TEET, DIL NEW, NEW 

note: NEW is the designation of the new data file. 

and the original DATA DECK is replaced with the one card: 


7 ADD, P NEW 
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r 


( PFIN 

( V FIN ~ (p 

( DATA DECK 


( PXQT DE 

© 

Optional Insert — ^ f LIB SYS$*MSFCFOR$ . ( 3 ) 

Pmap, inde.de 


f P PRT, T 


( PFREE NOZDS 

© 

f PCOP1N NOZDS, TPF$. 

© 

( PASG, T TPF$. f F///800 

® 

( PFREETPF$. 

_ ® 

( Prewind nozds 

® 

f PASG, TF NOZDS, T, 1177 5 

@ 


PRUN CONTROL, CARD — This is the first card in the run stream. 


(D 

0 ) 

0 ) 

© 

© 

© 

© 

© 


The analysis is finished 
Program is executed 

System subroutines are made available for program use 
The program is mapped 

A listing of subroutine names in program is printed 
The program tape is removed from tape reader 
Copies the program tape into computer storage 
Assigns sufficient computer storage for the program tape 
Returns the tape to its beginning 
Assigns the program tape 


CONTROL CARD SET-UP 


% 
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Once the new program tape is made containing the two data files, the 
program user has the option of using either of the two data files. This pro- 
cess can be continued until several data files are available for the program 
user. The data file chosen will be a function of the nozzle design problem 
being considered and can be called by replacing the DATA DECK with the 
one card 


^ PADD, P OPTION 

where OPTION is the name of the data file chosen. 

4.3 OUTPUT FOR' 'AT 

This section describes the printed output as well as the plotted output 
for the computer code. The example output presented is for the model nozzle 
design analysis performed for the wind tunnel test IA-604 planned for the 
NASA-Ames Research Center 11 x 11 wind tunnel (Ref. 1-1) and is the result 
of the input data presented in Table 4-2. 

The computer code output is organized so that the initial page contains 
a description of the model nozzle design problem being considered. The 
problem description is presented in such a manner that most of the input data 
are listed. If N PRINT has been set equal to 1, the second and subsequent pages 
will contain a listing of the intermediate data showing the results of the calcu- 
lations performed in Section 3.3, Steps 12 and 14 (Fig. 3-7). A typical print- 
out for each of the two steps are presented to demonst rate the output for each 
step. The final pages of the computer code output contain the final "envelope" 
of model nozzle designs, and corresponding power sweep operating character- 
istics of representative model nozzles for each of the desired fligh* - mcFtions 
to be simulated in the wind tunnel test. Numbered flags on the exariuo* print- 
out sheets correspond to the numbered comments in the following description 
of the printout. The calculations are performed in English units. 
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GROUP \ - MODEL NOZZLE DESIGN PROBLEM ’'SCRIPTION 


(T) This corresponds to the value of NMACH input on Card 1, 
and is equal to the anticipated number of scheduled flight 
Mach numbers to be simulated during testing. 

(2) Identifies the form of the similarity parameter used in 
the nozzle design analysis. 


Items (5) through (7) identify the prototype p- . ne similarity data for the 
rocket engine that the model nozzle is to simulate. Figure 3-2 presents a 
pictorial definition of these parameters. The values of M^, 6. and y . are the 
plume boundary properties corresponding to each value of P /P^. These 
properties are determined by means of a numerical integration through a 
Prandtl- Meye r expansion. 


(5) PC/PB: Prototype pressure ratio along the plume boundary 

(P /P, ) where P is the nozzle chamber pressure and P. is 
c' b c r b 

the vehicle base pressure. 


(?) MJET: Prototype plume boundary Mach number (M^). 

(§) DELJ: Prototype plume boundary initial expansion angle 


(6) GAMAJ: Prototype ratio of specific heats on the plume 
boundary (V ) . 

( 7 ) SP: This is the value of similarity parameter corresponding 

to each value of P /P, that results when using the values of 
c b 

A and B specified in item © . 


Items (§) through (fi) identify the flight conditions that the model nozzle will 
be designed to simulate during the wind tunnel tests. 


(8) MINF: Scheduled flight Mach number to be simulated during 
the wind tunnel test (M^), The values of P^/P^ , P c and P^ 

are those flight conditions corresponding to each value of M 
in the ascent trajectory. 
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PB/PI: Predicted value of Lise pressure ratio (P^/P ). 

PCs Prototype nozzle chamber pressure (P c ,psia). 

Pis Prototype freestream static pressure (P , psia). 

SPNOM: Value of similarity parameter required to ensure 
base pressure matching corresponding to each value of M 


GROUP 2 - OPTIONAL NOZZLE "FAMILY" DATA (not printed if NPRINT = 0) 


Items (1^3) through (23) define a "family" of candidate model nozzles which 
have been designed to ensure base pressure matching for " one" specified 
set of flight conditions . Each candidate nozzle in this family is defined by 
t> e and M and is designed based on a constant value of P £ (usually set equal 
to 50% of the maximum available air pressure supply). The nozzle flow 
separation constraint of Section 3.3, Step 1, is the only constraint considered 
at this point in the analysis. 


One "family" of candidate model nozzles is presented for each scheduled 
flight Mach number to be simulated during testing. 


(U) WIND TUNNEL PI: Wind tunnel free stream static 

pressure (psia). This is not necessarily equal to the 
prototype freestream static pressure. 



LIP ANGLE: Model nozzle exit plane wall angle 
deg) 



MEXIT: Model nozzle exit plane Mach number (M ^ 
dimensionless) c 


(l6) PC/PI: Ratio of model nozzle chamber pressure to 
wind tunnel static pressure. (This ratio is constant 
for a specified set of flight conditions.) 



PE/PB: Ratio of model nozzle exit plane static pres- 
sure to vehicle base region pressure. Nozzles with 
a P / P^ less than 0.60 will experience nozzle flow 

separation and are not considered further in the analysis. 


4-18 


LOCKHEED • HUNTSVILLE RESEARCH & ENGINEERING CENTER 


LMSC-HREC TR D784 ill 


(Th) PC/PBj Ratio o£ model nozzle chamber pressure to vehicle 
base region pressure. (This ratio is constant for a specified 
set of flight conditions.) 



PC/PE: Ratio of model nozzle chamber pressure to exit 
plane static pressure. 


(20) MJET: Plume boundary Mach number determined as a func- 
tion of PC/PB (Mj, dimensionless). 

(zT) AE/'A*: Model nozzle area ratio (A e /A*, dimensionless). 



DELJ: Initial plume expansion angle measured from th_ 
model nozzle centerline (6., deg). 



Lowest value of exit plane Macli number in the iterative 
solution that results in a PE/PB less than 0.60. Model 
nozzles with M >_ (23) are not considered further in 

the analysis. 


C.ROUP 3 - OPTIONAL NOZZLE "ENVELOPE" DATA (not printed if 
N PRINT - 0; not plotted if NPLOT - 0). 

The data presented in Group 3 define an "envelope" of candidate model nozzles 
that will: (1) ensure base pressure matching for " one " specified set of flight 
conditions to be simulated, and (2) meet all of the constraints that must be 
considered in the nozzle design study for that set of flight conditions. Each 
"family" of nozzles presented in the Group 3 data are those nozzles in each 
"family" of the Group 2 data that satisfy the model nozzle design constraint 


SP . . 5 .0 < SP . . SP - 5.0 

minimum — nominal — maximum 


where SP . - SPNOM corresponding to each flight Mach number, 

nominal c b ° 

SP . . corresponds to that value of P that will result in flow separation 

minimum r c r 

in the model nozzle, and SP corresponds to either the maximum 

maximum r 
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auxiliary air pressure supply or maximum mass flow rate capability of the 
test facility . 


One "envelope" of candidate model nozzles is presented for each scheduled 
flight Mach number to be simulated during testing. 


Figure 4-2 is a graphical representation of the Group 3 data. These plots 
are made automatically if NPLOT - 1 . 


(£ 4 ) SPMIN: Minimum obtainable value of similarity parameter 
for the given model nozzle and specified flight conditions. 
Corresponds to that value of P (PCMIN) that will result in 
flow separation in the model nozzle. 

SPMAX: Maximum obtainable value of similarity parameter 
for the given model nozzle and specified flight conditions. 
Corresponds to PCMAX . 

(Zb) PCMAX: Maximum allowable value of model nozzle chamber 
pressure. PCMAX is set equal to the lesser of either the 
maximum auxiliary air pressure supply or the chamber 
pressure corresponding to the maximum mass flow rate 
capability of the test facility . 


GROUP 4 - FINAL MODEL NOZZLE ENVELOPE DEFINITION AND RANGE 
OF APPLICATION 


Items (27) and (28) are the coordinates which define the "final envelope" 
of model nozzles that may be used in the test program. Each model nozzle 
in this envelope will: (1 ) ensure base pressure matching for "all" specified 
sets of flight conditions to be simulated, and ( 2 ) meet all of the constraints 
that must be considered in the nozzle design study. Figure 4-1 presents a 
computer code plot of the final envelope of model nozzles that may be used 
in wind tunnel test LA-604. A plot of the final envelope of model nozzles is 
made automatically independent of the value of NPLOT. 
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(29) Schedule of flight Mach numbers that may be simulated 

with the "final envelope" of model nozzles defined above - 
The power sweep operating characteristics will be pre- 
sented in the Group 5 data for each of these scheduled 
flight Mach numbers . 


GROUP 5 - POWER SWEEP OPERATING CHARACTERISTICS 

Items (So) through (^ 5 ) (fi) through (li) , (8) , (£ 4 ) and (£ 5 ) define the power 
sweep operating characteristics of five representative model nozzles that 
exist in the "final envelope" of Fig. 4-1 for each of the scheduled flight Mach 
numbers to be simulated. 


® 

@ 

@ 

© 


NOZZLE: Nozzle identification. See Fig. 3-8 for the 
relative location of eftch nozzle in the "final envelope." 

MEXIT: Model nozzle exit plane Mach number (M e » 
dimensionless ) . 

LIP ANGLE: Model nozzle exit plane wall angle 
(0^, deg) . 

MINE: Scheduled flight Mach number to be simulated 
during the wind tunnel test (M^ ). 

SPMIN: Minimum obtainable value of similarity param- 
eter for the given model nozzle and specified flight 
conditions . 

PCM1N: Value of model nozzle chamber pressure 
corresponding to SPMIN (psia). 

SPNOM: Value of similarity parameter required to 
ensure base pressure matching for the given model 
nozzle and specified flight conditions. 

PCNOM: Value of model nozzle chamber pressure 
corresponding to SPNOM (psia). 

SPMAX: Maximum obtainable value of similarity 
parameter for the given model nozzle and specified 
flight conditions. 

PCMAX: Value of model nozzle chamber pressure 
corresponding to SPMAX (psia). 
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(Tl) WIND TUNNEL PINF: Wind tunnel freestre&m static 
pressure (psia). 

(34) DSTAR: Model nozzle throat diameter (in.) 

(35) DEX1T: Model nozzle exit diameter (in.). 
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Example Computer Code Printout and Plots 
for the Model No/.zle Design Analysis Per- 
formed for the Wind Tunnel Test IA-604 
Planned for the NASA/Ames Research Center 
11 x 11 Wind Tunnel (Ref. 1-1) 
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Appendix 

A LISTING OF A COMPUTER CODE USED 
TO DESIGN MODEL NOZZLES WHICH 
MEET MSEC BASE PRESSURE 
SIMILARITY PA >. A METER CRITERIA 
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SUBROUTINE ARATIOIGAMWA ,KM,AR) 
G«|:CtMI>|.0 
GP1:GAMMA«1,0 
CKP:GP1/I2.0#GH1 * 

AT2.0/GP1 

6:1.0*IGM1/2.0>pxmpxh 
AR: | IABBIppCXP|/Xm 
RFTURN 
CM) 


SUBBOUT INI CHICK IXD,NN07,NP,XL,XR,KI 
DIMENSION KOI 20, 25,201, NPI20.2S) 

xb:u. 

XLJ 1 OOOOOOOOOCOO .0 
00 107 Jil ,NN0Z 
N:NP!l ,M I 

1MN.IO.OIGO TO 102 
00 103 j:l,N 

1 r « XDIl t J»K I -*B I 100,100,101 

101 KB: XO 1 1 • J ,K I 

100 If IXL-XOI 1 , J,M I 103, 103 , 104 
104 XL=XOt 1 , J«K I 
1C 3 CONT INul 

102 CONTI NUE 
R 1 1 URN 

t NO 


SUB ROl'T INI C0MB01IJMAX) 

C0HM0N/CQ1/GAHAJ,GAHAE,GAWAI,PATM,A ,B 
CCMM0N/C02/PCMAX , PC PBMN , PC 
C0HM0N/CQ3/NMACH ,NSP , NN 0? , NPR IN T ,NPL 0 T 
C0MM0N/CO4/SPll2O>,PeOPTl?0),XHIt?0>,Pll2O) 

C0MMCN/CO5 /PToPl I20»,PCOPII20I,PCOPBI?0»,XHJI701 .PHJI20I 
COM HON/C Ofc/PCCPt 120,251 ,XHf 1 20, 251, PHI! 70, 2 5), AO AS 1 20, 251 . 

• OIL TAJ! 20,2 5 I .TMITL I 20, 251 , PC OP 1 170,25) , PCOPB 1 20 , 2 5 » 

1 f 0RMATI8E10.5) 

2 r ORMATI 1H0, 1 IN I 1 SUB WNC 8110.5) 

3 I ORHATl IX, 10110.51 

4 f OR H A T ! I 5 ) 

5 f ORHATl 2X, •HIM : ',110. 5,* WIND TUNNEL PI = *,E10.5,* PB/PI : 

• • , C 1C . 5 , ' 5PN0H r *,110.51 

b FORMAT! 1H0, IX, *MEXIT IS SUBSONIC AND TEMPORARILY SIT : 1.10*1 
7 FORMATI 1M0,1 Xi*L 1P ANGLE H|XIT PC/PI Pt/PB RC/PB 

S fC/PC MJLT At /A* OELJ *) 

6 FORMAT! INC, IX, ‘FLOW SEPARATION OCCURS FOR HCXIT GREATER THA" *, 

• 110.5,' WHC N IINF - *,C 10.51 

9 FORMAT! IhI, IX, *TH| FOLLOWING IS A N027LC FAMILY OERlVCO TO SIMULAT 
•t THE FLIGHT CONDITIONS') 

JMAX=25 

00 130 1:1 , NMACH 

CP FREE STREAM PRESSURE RATIO 

CALL PR A T 10 I GAMA I, XMI II I.PTOPIII )) 

CP FREE STREAM STATIC PRESSURE 

PHI |:PATM/PTOPI II ) 

IFINPRINT.EO.OIGO TO 15 
WR1TEI6,?) 

WRITE 16 ,5)XMI II > ,PI II I.PB0P1 II »,SP1 II ) 

WR I T 1 1 6 , 7 I 

15 CONTINUE 
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Subroutine COMBO 3 (Cont'd) 

PCOPIII ISPC/PIII I 
PCOPBII ISPCOPIII I/P80P1 m 
C* MACH NUMBER ON the jet boundary 
CALL MACMtGAHAJ.PCOPBtl I.XMJUI I 
C* JET BOUNDARY ISENTROPIC EXPANSION ANGLE 
CALL PRNOIGAMAJ.XMJIII.PMJIIII 
t* 

C* THE FOLLOWING STATEMENTS CALCULATE NOZZLE LIP ANGLE FOR 
C* A DESIGNATED PC/PE 

C* 

PCOP£tI,lir$.f 
00 110 JSltJMAX 

C* MACH NUMBER ON THE NOZZLE EXIT PLANE 
CALI MACM46AM £ tPCOPEU tJI »XME( I»JI I 
CP NOZZLE LIP 1SCNTROPIC EXPANSION ANGLE 

CALL PRNOtGAMAE. XME 1 1 , J I , PNC i I • J 1 1 
C* NOZZLE EXIT PLANE AREA RATIO AE/A* 

CALL ARATIO IGAMAE ,XMEII , J I , AOAS 11 , J I I 
C* PLUME BOUNOARV INITIAL EXPANSION ANGLE 

OELTAJI I,JI = ISP1 1 1 l/XMJt 1 1 I PI XME 1 1 ( Jl PPAeGAMAjeeB I 
C* NOZZLE LIP ANGLE 

THCTLII •J»:0ELTAJI1,J|«PNEI I tJl’PMJH » 

PCOPIII .J1=PB0PI II lePCOPBU l/PCOPEtl, Jl 

PEOPBII , JISPC0P8 III/PCOPE II .Jl 
IFINPRINT.EQ.OIGO TO 20 

WRITEI6, JITHE’LI I. J) .XMEII, JI.PCOPI 1 1 I , PEOPB 1 1 , J I , 
SPC0PBI1 I.PCOPEII .JI.XMJIl I.AOASll.JI.OELTAJtl.JI 
20 CONTINUE 

IFIPEOPBtl.JI .LE .PEPBMNIGO TO 120 
IFIXMU l.JI .GT.3.5IGO TO 30 
PCOPEII . J«lliPC0PEII,JI*S.0 
CO TO 100 
30 CONTINUE 

PCOPf (l.J«l):PC0PC!l.JI«20.0 
ICO CONTINUE 
11C CONTINUE 
GO TO 1ZS 
120 CONTINUE 
JPAX r J 
125 CONTINUE 

IFINPRINT.EQ.CIGO TO 130 
URITEI6.BIXMEII.JI.XHK1I 
130 CONTINUE 
RETURN 
END 


SUBROUTINE OESIGN 
CALL INPUT 
CALL COMB03IJHAX) 
CALL SPRNGE I JHAX I 
CALL FINAL 
RETURN 
£ NO 
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SUBROUT I HE FINAL 

C* THIS SUBROUTINE DETERMINES THE FINAL ENVELOPE OF NOZZLES THAT HAY 

C* IN THE TEST PROGRAM SUBJECT TO HASS FLOM ANO FLOW SEPARATION CONST 

COHMON/COt YGAMA J ,BAMAE ■ 6AHA 1 ,PA TH, A , B 
COMMON/ CO J/NHACH.NSP, NNOZ, NPRINT.NPLOT 
COHNON/CO«/SPH?OI ,PBOPI I ZO I, XM 1 1 20 I .PI! 20 1 
C0HH0N/CO«/XHCG!2O,2S,2OI,TMCTLOI?0,25,2OI ( NPI2O,?&> 
C0MHON/CO9/HCDRPII2,20l ,XLINC(12I«XF1NAII12I,VF1N*LI 121 
COHHON/TEST l/XHCL,XMER, THLB iTHLT ,THUY ,THUB,HHAX 
OIHENSION ThCTLII20I,ThCTLUI20> 

0 1HENSI ON CNVCL! A0l,vTlT||2l,X1I1ll2l 

OATA IENVCL!II,1:|,BDI/'FINAL envelope of hooel nozzles THAT HAY b 
•E USED IN THE TEST PROGRAM WHICH HECT THE SIHILARITV PARANlUR 
SCR1TCR1A, NASS FLOW ANO FLOW SEPARATION RESTRICTIONS IS PLOTTED B 
SELOW. THE SIHILARITV PARAHETER USEO IN THE NOZZLE DESIGN STUDY 
SIS OCFINEO BY THE REL AT lONSHIPSP^M JE T *DC L J/ IMC X I T**A*6 AHA J#*B 1 . WH 
SERE A: ANO B - . THE HOOEL NOZZLES REPRESENTED BEL 

SOW HAVRE USED TO S1HULATE THE FLIGHT CONDITIONS THAT EXIST FOR AL 
SL FREE STREAM MACH NUMBERS TO BE TESTED. */ 

DATA BLANH/bH / 

1 FORMAT! |H1. IX, 'THE FINAL ENVELOPE OF HOOEL NOZZLES THAT HAY BE USE 
SO IN THE TEST PROGRAM WHICH NEE T THE SINIL ARI T Y • tit , *P ARAHE TER CR1 
ATCRIA, HASS FLOW AND FLOW SEPARATION RESTRICTIONS ARE DEFINED BY T 
SHE COORDINATES' I 

2 FORMAT! 1H0.13X, 'NOZZLE CUT MACH NUMBER ', 2** , 'NOZZLE LIP ANGLE (DC 
• Gl'l 

I FORMAT! IX, l9X,E10.S,SbX,E10. 51 

k FORMAT! IMO, IX, 'THE FLIGHT MACH NUMBERS THAT MAY BT SIMULATED WITH 
STH1S ENVELOPE OF NOZZLES ARC • I 
S FORHAT!?X,'FR[£ STREAM MAfM NUMBERS :',10C10.S> 

C* THE FOLLOWING STATEMENTS DETERMINE THE UPPER LIMIT ON NOZZLE EXIT 
xmcr :iocioro.o 

XMEL=r.O 

DO 10? 1=1, NNOZ 

NrNPlI.ll 

lFiXMCR-XMCOI 1,N,1 I I 102 ,102 , 101 

101 XMER:XMCOI I ,N,1 I 

102 CONTINUE 

C« THE FOLLOWING STATEMENTS DETERMINE THE LOWER LIMIT ON NOZZLE EXIT 
UMAX :NM ACH 
10S CONTINUE 

DO IDS 1=1, NNOZ 

IF 1 XMrO I 1, 1 ,HMAX l-XMEL 1 105 , 105, 10* 

1TN XMCL^XMEOI 1 , 1 .UMAX I 

105 CONTINUE 

IF! XMfL.LC.XMCRIGO TO ICb 
KMAX rXMAX-l 
GO TO 103 

l* THE FOLLOWING STATEMENTS DETERMINE THE UPPER ANO LOWER LIMITS OT N 

C* LIP ANGLE CORRESPONDING TO THE LOWER LIMIT AND UPPER LIMIT OF NOZZ 

C* MACH NUMBER RESPECTIVELY 
10b CONT 1NUC 
JJ:0 

00 110 1:1, NNOZ 
N:NP 1 1 , K MAX I 
00 109 J: 1 ,N 

IFIABSI XMEOI I «J,HHAXI*XM£L I ,LC. • 0 1 1 GO TO 10B 
GO to 109 

106 CONT 1NUC 
jjrjj. i 

TMETLL! JJ!:TMETL0I1,J,KMc* I 
GO TO I 10 

109 CONTINUE 

110 CONTINUE 
J JM A X : J J 
KM:0 
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OZZLE 
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Subroutine FINAL (Cont'd) 

00 IIS 1 — 1 » NN02 
NiNPl 1,11 
00 114 JU.N 

IF! 6BS I INC 01 1 » J » 1 1 -MUCH I «IC • *01 160 10 IIS 
GO 10 111 
US CONTINUE 
kr:rr*1 

THE TlUIRRIUMCTLOt l.J.l I 
GO TO 11$ 
lit CONTINUE 
11$ CONTINUE 
rrn*x:r : 

ThlBUOO' JO.O 
TMLTrC.C 

thug : iooooo *o 

iHUTrc.O 

00 12C I :i . JJNAX 
IF I THCTLLI I l-THL Tl 11 7, 1 17, 1 16 
lit Tml T:»mCTLL 1 1 » 

117 IF* T Hi B -THCTLLI I II 119,1 19.116 

118 T HLB : IHE TLL I I I 

119 C ON 7 I NUC 

120 CONTINUE 

00 12* UI.RKNAX 
IFITMEUUtl I-THUTII22.122.I21 

121 T HU T : THE TLU I I 1 

122 IFtTMUB-THCTLUlI 11124,124,123 

123 T HUB : THE TLU 1 1 I 
12<* CONTINUE 

12$ CONTINUE 
NRITEtb.il 
MRITEI6.2I 
NRITCI6.3IXNCL.THLB 
NRITCI6.JIXNEL.THLT 
NRlTCIb.3IXHCR.THUT 
NRITEIb.3IXHER.THUB 
NRI TC lb .4 I 

NRITCI6.5IIXM till, 1:1, RHAXI 
CNCOOE I 696.CNVCL 1531 1* 

CNCOCEI 696.CNVCLIS6I IB 
696 F0RHAT(F6»3I 

CALL CHSI2VI2.2I 

CALL RI TE2V ISO, 1000.102 3.90.2.460.1 .CNVCL. ICR I 
C* LOAO LABELS FOR AXES 
DO 699 11:1.12 
VT1TIII | :YF INAL I II I 

699 XTIT I 11 UXFINALI III 
00 700 KHM ,12 
JR:RR 

IF I XT1T I JR I .EC. BLANK IGO TO 701 

700 CONTINUE 

701 CONTINUE 
jr: JK-2 

1 XSHF =5 1 2- JR*54 
L* St? UP GRID ORIGIN 

CALL SClHlV I 30, 30,34,2501 

lNt:0 

INX:0 

THIN: THLB 

YNAX :THUT 

"HI N SXHEL 

XHAX :XHER 
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Subroutine FINAL (Gont'd) 


CALL XYLlMI YM1N, YMAX, IV Y.AMIN.AMAX.DY ,NV,NYL tNY2l 

VMlN:AHIN 

yxiciNii 

CALL XYLiHtXNIN.XNAX , IUX,AHIN,AMAX,DX »NX *NXL »NXZ I 

ih]n:amin 

XMAXIAMAX 

SCI UP GRID 
CALL CHSIZVI2 f 2l 

CALL r»RIDI«<2(XMINtXNAX ,VMIN , VRAX tDX ,0V tNX ,NV ,NXL ,NYL ,NX2 ,NYZ I 
PRINT VERTICAL SCALE LABEL 

CALL R1TE2VU 2,200, 102J , I AO , 2 ,A8 , 1 , V T IT , IER I 
PRINT HORIZONTAL SCALE LABEL 

CALL RITC2VI 1 XSHF , lb , 1023,90 , 2, AS , 1 ,X TIT, ICR I 
PLOT FINAL NOZZLE ENVELOPE 

ixhcl:nxvixheli 

ixhcr:nxvixhepi 

1 THL TiNVVITMLTI 
ITHLB-NVVITHLbl 
I THUT :N Y V I THUT l 
ITHUB:NVV<THUI’I 

CALL LINEVt IXHEl .ITMLB, TXHEL.ITHLTI 
CALL lINCVtIXHCL , I THL T , 1 XME R , I THUT » 

CALL L1NCVI IX**CR l ITHUTtIXPER|ITHUB> 

CALL LINCV(1X M CR,I THUB «IXMEL»ITHLB* 

CALL TEST 

RETURN 

END 


SUBROUTINE FIRST tNNOZI 
EXTERNAL TABLAV, TABLIV 
COHHON/XIN/TI TLE 181 
COMMON/ CD/N0Z0«10 t 20 » 

DIMENSION HE AD 1 6 I 

0A1A ME AD/6HSVHB0L ,bH N072 ,6HLE FAM.bHILY 1 0 , 6HE N T I F ! , bH C AT ION / 
CALL F R AME V 1 0 l 
CALL CH S12V 1 1 « 3 I 

CALL RI TC2V« ISO, 995, 102 3 , 90 , 2, AS, 1. T1 TLE, ICR» 

PRINT TABLE HEADING 
CALL CHSIZVI2.2) 

CALL R! TE2V»lAA,b0O,1023,9O,2,3b,l,HCAO,IERI 

SYMBOL TABLE 
I X: 1 7 A 
I Y-SfcB 
NY: IY-10 
I JMI 1 

CALL CHSI2VI2.21 
DC IS 1:1 »NNO 7 

CALL VCHARVt90,2,IX,NY,IJK,TABLAV» 

CALL RITSTVI 13,12,TABL1V» 

NA:N V*b 

CALL PITC2Vt290,NA,lO23,9O,2,bO,l,NOZD( 1 ,1 » , IERT 
I y: I Y -20 
NYHY-10 
IS I JM: 1 JM* 1 
RETURN 
END 
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SUBBOUT INI GRAPH IYMAX,YMJN,XNAX,XN1N,XC,0CN,XT IT ,VT IT , NP.NNOI ,K I 

C 

C* • SUBBOUT INC GRAPH SITS UP THE GRIDS AND CONTROLS PLOTTING 
C 

EXTERNAL TABLIV 
COHMON/XHS/IXS 
CCMMON/XIH/TI TLt 181 

C 0MH0N/C04/SP 1 1 SOI . POOP II SO I, XH II SOI. PI I SO I 
696 FORMAT I F6, 3 1 

DIMENSION XT1TI12I.YTIT US I .SYMBOL I SSI 

OA7A ISVMBOLIII.IbI.SSIF'ThC SYMBOLS PLOTTEO BELOW REPRESENT MODEL 
S NOSSLES THATMAV BE USED TO SINUL»TE THE FLIGHT CONDITIONS THAT 
S EXIST FOB A FREE STREAM MACH NUMBERS V 

DIMENSION NPISO.2Sl.XDISO.SSI.DCNI20.SSI 
C 

C** ADVANCE FRAME 

CALLFRAMCVIOI 

C 

C*« SETUP GRID ORIGIN 

CALL SETMIVI3 ,30,68,1341 
YSAVVSVMAX 
ysava: ymin 
X SA V V "X max 
xsava:xmin 
c 

Iwv:q 

lwx:0 

CALLXYLIMIVMIN.VHAX.IWV , AM I N , AH AX ,0Y , NY ,NYL ,NY2 I 
YMIN SAM IN 
VMAX SAMAX 

CALLXYLIMIXMIN.XMAX.IWX.AMIN.AHAX.OX.NX.NXL.NXSI 

XMINSAM1N 

XMAXSAMAX 

c 

C** SET up grio 

CALL CHSI2V 4 2 ,2 I 

CALL RITSTVI 1S.S6.TABL1VI 

CALL GRIOIVIS.XMIN.XMAX , YMl N , YMA X ,0 X , DY .NX , NY ,NXL ,NYL .NX2.NVSI 
YMAX SYS A VV 
YMINSVSAVA 
XMAXSXSAVV 
XMINSXSAVA 
C 

c 

C*« PRINT VERTICAL LABEL 

CALL RITC2V(12,404,1023,IB0,2,36,I,VTIT,1ER> 

C*« PRINT HORIZONTAL SCALE TITLE 

C ALL RITE2V4 1X5,50,1023,90,2 ,36,1 ,XT IT .Itfil 
C 

C** PRINT FIGURE TITLE 
CALL CHS IZV I 3 ,3 1 

CALL R1 TE2V 1 1 50. 10. 1023,90 ,2 .OB. 1. TITLE, IERI 
C* * PRINT DEFINITIVE NOTE 
CALL CHSIZVIS.2I 
CALL R1TSTVI16.26.TABL1VI 
ENCODE I 696. SYMBOL I 25 1 IXMItKI 

CALL RI TE2VI5D, 1000, 102 3, 90, 2,150,1, SYMBOL, IERI 
43 CONTINUE 

I JMr 1 

DO 299 INK: 1 ,NN02 
NPJ:NP| INK, HI 

00 301 Jj:i ,NPJ 

C ALLPOINT XI XOIINK, JJI ,OCN(INK,JJI ,1 JM I 
301 CONTINUE 

1 JM : I JM ♦ 1 
299 CONTINUE 

RETURN 

END 
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Data File HPM 

* I I 1 5 

.02 SCALE MPM N022LE CUIT HACK NUMBER 
.02 SCALE MPM N022LC LIP ANGLE 

.02 SCALE HPM N022LC OESIGN ANALYSIS TEST IA-609 
N022LC LIP ANGLE 


17 

37.69 

7.6599 

10.990 

1.2591 

97.97 

2.7869 

19.258 

1.2685 

61.99 

2.9223 

18.025 

1.2805 

80.99 

3. 0761 

21.793 

1.2893 

107.32 

3.2393 

25.561 

1.2887 

199. 28 

3.9131 

29.329 

1.2939 

196.99 

3.600S 

33.096 

1.2997 

273.55 

3.8026 

36.669 

1.3062 

387.20 

9.021? 

90.632 

1.3129 

559.23 

9.2585 

99. 395 

1.3189 

876. 99 

9.5209 

98.167 

1.3298 

1755.1 

9.8153 

51.935 

1.3326 

1960.9 

5.1999 

55.703 

1.3906 

3167.9 

5.532C 

59.971 

1.3985 

5317.3 

5.96 76 

63.237 

1.3559 

9336.8 

6.9766 

67.006 

1.3627 

17781.0 

7.0879 

70.779 

1.3687 

15U0.0 

750.0 

0.6C 

19.70 

.0? 

199.6« 

27.' 

500.0 

.597 

.89 

72C.O 

10.592 


1 .AO 


N027LC FAMILY DERIVED FOR FLIGHT MACH NO. : 0.597 

N027LC EXIT MACH NUMBER MINF=0.59T 

.796 .825 651.9 B.991 

N022LC FAMILY OCRIVCO FOR FLIGHT MACH NO. S 0.796 

N072LC EXIT MACH NUMBER MINF:0.796 

.90 .77 629.5 7.925 

N077LE FAMILY OERIVCO FOR FLIGHT MACH NO. : 0.90 


1 .90 


MINFS0.90 
6.8 77 

FLIGHT MACH NO. 
MINF:0. 95 
5.725 


= 0.95 


: 1.098 


NC22LC EXIT MACH NUMBER 
.95 .675 610.2 

N022LC FAMILY DERIVED FOR 
N022LE EXIT MACH NUMBER 
1.096 .61 580. 0 

N022LC FAMILY DERIVCO FOR FLIGHT MACH NO. 

N072LE EXIT MACH NUMBER HINF=1.098 

1.10 .69 571.3 5.977 

N077LE FAMILY DERIVED FOR FLIGHT MACH NO. = 1.10 

N022LE EXIT MACH NUMBER HINFrl.10 

1.196 .65 570.0 9.877 

N022LC FAMILY DERIVED FOR FLIGHT MACH NO. = 1.196 

N077LC EXIT MACH NUMBER M!NF:1.198 

1.299 .69 57b. 2 9.105 

N027LE FAMILY DERIVED FCR FLIGHT MACH NO. S 1.299 

N027LE EXIT MACH NUMBER MINF=1.299 

1.903 .77 586.0 3.276 

N077LE FAMILY OCRIVCO FOR FLIGHT MACH NO. : 1.903 

N077LE EXIT MACH NUMBER MINF=1.903 


1.90 
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10 


C* 

c* 

c* 

c* 

c* 

t* 

c* 

c* 

c* 

c* 

c* 

c* 


SUBROUTINE INPUT 

C 0 MH 0 N/C 01 /GAHAJ .BAMAC , GAMA I ,PA TM , A ,0 
C 0 MM 0 N/C 02 /PCMA!) .PCPBMN.PC 
C 0 MM 0 N/C 03 /NHACM,NSP.NN 0 I,NPRINT .NPLOT 

common/con/sp i 1 20 i ,pbopii 2 oi,iph 2 oi,ph 20 i 

C 0 MMON/CO 7 /SCALC «PC II TO .WMAI.TO 

CuMM 0 N/e 0 */HC 0 RPI 12 . 20 l .III NCI 121 .IF INAL I I 2 I.VFINALI I? I 
COMMON/ T CST 2 /NN 07 F 
COMMON/ II M/ TITLE IBI 
C 0 MM 0 N/C 0 /N 070 I 10.201 

OIMCNSION PPCOPB I JO) .PXMJI 30 ! .PDCLJI SOI .PG 5 MAJI SO » ,PSP I JO » ,PPC I 201 
t, PPM 20 I 

1 TOR MAT MIS | 

2 F ORMA T I BE 10 • S I 

3 FORMAT) IHl, IX, * 4 C 6 ;,N A NCR N 022 LC DESIGN STUOY. IT IS ANTICIPATED 

• THAT INC NOZZLE WILL BE DESIGNED FOR M 2 .’ FREE STREAM MACH NUMBC 

• RS. 1 ! 

13 rORMATI 1 HQ. 1 I.MhC SIMILARITY PARAMETER USED IN THIS N 022 LC 0 CSI 6 N 

• STUDY IS DEFINED BY THE PEL A Tl ONSM I P •/ 21 , * SP S M J£ T *DEL J / 1 ME X I T •• 
SA*GANAJ»» 6 I . KMERC A : •,F 10 . 5 ,> ANO B S '.CIO.S! 

A FORMAT I 1 H 0 , SI. 'THE N 022 M U ILL BE DESIGNED BASCO UPON THC'/ 6 l,'F 0 l 
SLOWING PROTOTYPE PLUME SIMILARITY QATAM 
S FORMATIII.IOEIO.SI 
b rORMATI 1 ND. 1 I,> PC/PB M JC T 

7 FORMAT! 1 H 0 . 6 X, 'THE N 022 LC WILL BE 
AOlLOwlNG FLIGHT CONDITIONS'! 
e FORMATI 1 H 0 . 1 X,' MINE PB/P 1 PC PI 

9 FORMAT) 10 A 6 ) 

12 FORMAT) 12 A 6 .Al. 12 ! 

PCMAX- MAXIMUM ALLOWABLE N 022 LE CMAM 6 C R PRESSURE IPSIA! 

PC:H 1 ORANGE VALUE OF N 022 LE CHAMBER PRCSSURC USUALLY : 
PIPBMSrCXPERlMENTALLY OPSERVEO VALUE OF PCXIT/PBASC AT 
SEPARATION OCCURS IN THE NOZZLE 
PATMITOTAL PRCSSURC OF WIND TUNNEL IPSIA! 

GAMA j:RAT 10 OF SPECIFIC HEATS OF THE N 072 LC EXHAUST AT 
i> AM A £ - RATIO OF SPECIFIC HEATS OF THE N 022 LC CIHAUST IN 
PLANE AT The N 027 LE LIP 
GAMA I :R A 1 10 OF SPECIFIC HEATS OF AIR 

NMACh:ANT 1 CIPATE 0 NUMBER OF FREE STREAM MACH NUMBERS 
THAT THE N 072 LC WILL BE DCSIGNCO FOR 

NPMR AY: NUMBER OF RAYS IN PROTOTYPE PR ANOT L -MEYER EXPANSION OATA 
READ IS, 1 INHACH.NSP.NPRINT.NPLOT , NN 07 F 
RE ADIS. 12 IXFINAL 
RLA 0 IS. 12 IVF 1 NAL 
IFlNPLOT.CO.OIGO TO 40 
READIS. 9 IITITLCI I I , I a 1 , R | 

REA 0 I 5 , 1 2 I X L 1 NE 
40 CONTINUE 
NNOZ^NMACH 
GO TO) 10 , 20 , 3 C 1 ,NSP 
CONTINUE 


OELJ 

DCSIGNCO 

PC 


gamaj SP 
SIMULATE '/IIX.MHE 


SPNOM 


PCMAX/ 2.0 
WHICH FLOW 


the 

thc 


PLUME B 
N 072 LC 


IN THE TEST P 


1 C 


2 D 


30 


35 


A:, 2 S 
B - 1 • 0 
GO TO 35 
CONTINUE 
A :. 2 S 
B :.5 

GO TO 35 
CONT INUC 
A : 0 . 0 
B:l.O 
CONTINUE 


OUNOARV 

EXIT 


rogham 
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Subroutine INPUT (Cont'd) 

*RITE(6,3INMACH 

WRITEU.ISIA.B 

REA0(5,2IPXHE 

REA0(S,1 INPRR* V 

WRITEI6,«I 

WRITE(6t6l 

DO 100 I:1,NPMRAV 

READ(5,2»PPC0PB( I I ,PXMJ ( 1 1 ,PQEL J( 1 1 ,(>6ANAJI1 1 
P&P I 1 ISPXMJI I l*P0CLJ( 1 I /(PXHC*«A«PGAMAJ( 1 l#*BI 
WRITE(6,S>PPC0PB(I>,PlHJ(ll,PCtLJ(I I , PG AMA J( 1 1 ,P$P I II 
PPCOPBI I »:ALOG(PPCOPB(I II 
100 CONTINUE 

RCADI5,2IPCMAX,PC,PEPBMN,PATN,GANAJ,6AMAE,GAHAI 

REA0IS,2ISCAL£,PEXIT0,VNAI,T0 

WRITE 16(71 

WRIT£(6,6I 

DO 200 ISI.NmACH 

RCAO IS, 2IXMI 11 1 tPBOPl (1 I.PPCt 1 |PP 1(11 
IFINPLOT.EQ.OIGO TO 110 
READIS ( 9IIN020(II,II,II=1,10I 
REAOI5,12l IHEORP (K , 1 1 # > 1,6 I 
110 CONTINUE 

PCPB:PPCm/IPBOPllII*PPIII 1 1 

PCPBsALOGIPCPBI 

00 250 j:i,NPMRAV 

JS:J 

IFIPPCOPBIJI.GT.PCPBIGO TO 260 
250 CONTINUE 
260 CONTINUE 
JTsJS-l 

FACTORS ( PCPB-PPC OPBI JTI 1/ 1 PPCOPBi JS » -PPCOPB 1 JT 1 1 
SPHI I:PSP(JT|*F ACT0R9IPSPI JSI-PSPI JT II 
WRITE <6, SIX HI (It .PBOPII II.PPCI1 I.PPIllI ,SP! (II 
200 CONTINUE 
RETURN 
END 


SUBROUTINE ITSUS (FOFV,V,SAVE,CONV,NTIMESI 
C THIS SUBROUTINE FROVIOES ITERATION CONTROL FOR ANT FUNCTION 

C SAV£(R'7I=ST0RAGE LOCATIONS FOR X AND FOFK 

C (CON VI 'CONVERGENCE CRITERIA 

C INTIMESI-HAX HUNGER OF ITERATIONS 

OIHENSIONSAVE 16 I 
NUSAVEI3I «.l 
F OF XCKXSAVE 1 8 I 

C F OF t ANO V ARE DUMMY INPUT ARGUMENTS 

FOFX cFOFV 
X:v 

C CHECK FOR CONVERGENCE 

IFIABS (FOFX l-CONV.LE.O. IGOTOIIO 
I T I HE - S A VE ( 1 I «. I 

C IT1ME CONTROLS THE TYPE CALCULATION TO BE PERFORMED 

C ITIMC:lt FIRST TIME THROUGH 

C ITIM£=2,P0S FIRST TIME THROUGH 

C ITI*E=J,NEG FIRST TIME THROUGH 

C ITIME=N, SOLUTION IS BRACKETED 

C ITIME=5, SOLUTION HAS CONVERGEO 

C ITIMES6, SOLUTION WILL NOT CONVERGE 
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Subroutine 1TSUB (Cont'd) 

6010(10, 30, 50, 701, ITIME 

initialize 
to N|S| 

ITINC:: 
rorxcKirorx 
save «•» srorxc* 

IMFOri.Ll.O. IC0 1050 
10 IFirori.Ll.O. 1601070 
I rtr or iCH.Gc.ror H601055 

SAVE (21 

X:X-2.*SAVCI2I 
601090 
55 SAVC(A|:x 
save ( 5i srorx 
1=X-SAVC(7I 
or ONE V API ABLE 

irorn-ruNCiicN which is driven io zero 

( II -VAR1A8LC WHICH IS ITERATIVELY SOLVED FOR 
(SAVE l-PROGRA* CONTROL 
SAvr«| I = IT1*'E 
S A VE I 2 I * X INC REHCN1 

SAVE I 5IXC0UN1CO DENOTING NlH ITERATION 
G01090 
50 I T 1 HE - 1 

IMF OF* .6T.0. IGOTOTO 
IF IFOFiCN.Li «FOF XI G0 105 5 
SAVE 121 r-i.*SAVC 12 I 
X1X«2.*SAVE 121 
G0T090 
55 SAVE Ifcl -* 

SAVE I7| rrOFX 
*:X*SAVC(2I 
601090 
70 I 1 1 HE XA 
Nl-SAVC I 31 

lnrori .L1.0. IGO T 0 75 
SAVE IN I IX 
SAVE « 5 I =FOr X 
GOT 060 
75 SAVCIFIXX 
SAVE (71 xFOFX 

RI Cl. NEW GUESS FOR X ACCORDING TO TYPE CALCULATION 
80 *XSAVC(4|-SAVE(5:*( ( S A V E ( 6 I -SA VC I A I I / (SAVE ( 7 I -SAVE (5111 
90 irtNl.GE.NTIHCSIGOIOlOO 
NI = N IM 
SAVE ! 3 I =NI 
GOT 0 120 
100 IUHE:6 
GOT 0 1 20 
I1C 1T1HEX5 
SAVE (AC* 

SAVE (51 IFOFB 
SAVE Ifcl =X 
SAVE 171 •tfofx 

120 SAVE ( 1 I XFLOAT (IT IHE I • .1 
Y :X 

RETURN 

END 
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SUGROUT INI MACH! GANNA, P, AMI 
GN106=»GANNA-I.I/SANNA 

GN|:gamna-1. 

INrSORT 12 •• IP**GM IOG*l • I /CM 1 1 
RETURN 
I NO 


MORE 

CONMON/XIN/TITLEIBI 

COHMON/CD/N020ll0,20t 

COMNON/COl/NMACH ,NSP,NN07,NPRINT ,NPLOT 
C0NN0N/C0t/XMCD«20,25,20»,TH£TL0I20,2S,20I,NP<20,25> 
CCNNON/C09/ME0R*«l2,2O> ,XLINE<12»,XF1NAL 1 1 2 I , VF I NAL 1 1 2 I 
COMMON/ IMS/ | XSMF 
DIMENSION XTITI12I,VT1T<12> 

INTEGER AOARY 122 • 

DIMENSION IAIISt20,2S> 

DATA BLANK/ AH / 

DATA IAOAPVI 1 1, !:|,S1/' I OS mm 1 HARO COPY REQUIRED •/ 

CALL IOENTI IOS.AOARVI 
CALL FRM0CEJ11 
00 19 1:1,20 
DO 18 j:l,I2 
HEORPt J,II:BLANK 
16 CONTINUE 
19 CONTINUE 
CALL DESIGN 

IT INPLOT .EO .OIGO TO 209 

C* LE AO FRAME IDENTIFIES N027LE OCSIGN PROBLEM AND N077LE FAMILY NUMB 

CALL r I RST I NN02 1 
00 2 UO K:i,NMACH 

DETCkMlNC PLOT LIMITS OF H0RI70NTAL AXIS 
CALL CHECKIXMEO, NN02 ,NP,XL,XR,MI 
C*« CALCULATE Y AXIS LIMITS 

CALL CMECKHHETLO,NNOZ,NP,Y8,YT,KI 
C** SET THE DEPENDENT VARIABLES 
DO 109 I : 1 , NN02 
NtNP I! , M I 
DO 10B j: 1 , N 
KAXISII ,J»:xmFD( I,J,K» 

1C8 CONTINUE 
109 CONTINUE 

C*« LOAO TITLE FOP ARES 
DO 699 11=1,1? 

YTITII1 irXLINUIIl 

699 XTIT (II 1 =ME DRPl I I, Ml 
DO TCP KM-'l ,1? 

JK : M H 

I F « XTIT I JMI .EO. BLANK 1G0 10 701 

700 CGNTINUC 

701 CONTINUE 
JK: JK -2 

IXSMf:512-JK*59 
C9* PLOT OATA 

CALL GRAPHIVT , YB , XR , XL , X AX I S , THE TLD 1 1 ,1 ,K I , XT I T , VT IT ,NP,NN07,K) 

200 CONTINUE 
209 CONTINUE 

CALL ENDJOB 

STOP 

END 
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E UNC T IONNXV 111 
CALLXSCLYHX.II, IERI 
NIV2|I 

• hubs 

C NO 


EUNCTIONNTEMTI 
CALL VSCLVI I T , IT, IERI 
NT V I It 
RETURN 
( NO 


SUBROUl I Nt POINTXI I.T.NSI 
t X T C RN A L I ABl *• V 
HSrIAPSINSI 
CALL CMSWVIJ.JI 
llUNXVt X I 
I trNtvi t ) 

NJ! : I x -6 
N» : I Y-9 

CALL VCMARWI9i 1 ,? t vx,Nt,nS.TA9LNVI 

CALL CHS 12 < - '.Si 

RETURN 

CSC 


SUBROUTINE PRATICIGAHHA ,XH,PR) 

OHIO 2: I GAPHA-l.T 1/2.0 

G0GN1 2GAHHA/ICANNA-1 .01 

PR:« 1 .0*GN102*XH*XH)**G0GNl 

RETURN 

END 


SUBROUTINE PRNDI GANNA, X*,PHI 
GPIGM1:IGAMNA«1. I/IGAHNA-1. I 
(i “ 1 G P 1 - 1 ./GP IGNI 
X*S1 :XH*XH-1 . 

E ArSCRT <GNIGP1*XMSI I 

TPiSORT I XHS1 I 

EC: SQR T IGP1GMII 

PP:rC*ATAN(EA I - A TANIF8 I 

PMrp**5T. 2957795 

RETURN 

END 
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SUKBOUT INC SPPNGI < JNXX I 

C'MON/COI/GAWAJ.BANAC.SANAI.PATN.A.B 

COHNON/C07/PCNAX .PIPBWN.PC 

C ON N ON / C 0 1 / NN AC H . N if , NN 0 1 . NP R I N T , NP L 0 1 

CONNON/COP/Sf 1 1 701 .PROP 1 1 70 I . XN 1 1 20 1 . PI I 20 1 

CONNON/ COS/P TOP I I20I.PC0PI I70I.PCOPBI20I .INJI20I.PNJI20I 

CONNON/ COA/PCOPC 120.251 .XNCI20,?SI.PNCI20,2SI.A0ASI20.25I. 

SOtlTAJI 20.29 I.THCTL 120. 79 I.PCOP I 170.251 .PCQPR 120. 25 1 
C ONN ON /COT/SC ALC .PCUTD.WNAX.T0 

C ONN ON/C OS /x N[ 0120. 25. 20I.THCTL0I 20 <25.201. NP 120. 251 
OINCNSION DC NON 1 70 .25 I . PC UOP 1 120 .25 I .fCUOfi l20.29l.INJNiK20.25l. 
SPnjnai I 20.751 *0C L JNX 120 .251 .if Nil 120 .251 .PCLOPBI 20 .75 1 . 
S1NJN|NI20.2SI.PnjNINI70.2SI .0CLJNNI20.25I.5PNINI20.25I 
I FORNATI BCIO.SI 
} FORNATI IX, IOC 10.51 

« fornatiihi.ix.’of tHC previously ocf i vco Nome faniucs. iMt rcu 

SOWING N0/7LLS NAY SC U5C0 TO S1NULATC THC FLUNT CONDITIONS*! 

& F 0 RNATI 2 X,*NINF : *,CIO.S.* NINO TUNNEL f| s *.CI0.5.* fil/fl : 
i • ,C 10.5.' Sf NON : * . C 10. 5 1 

6 FORNATI | HO I 

T F ORNAT I IHO. I*i* NC1IT LIP ANGLE Sf H|N SPNAX PCNAX *1 

CxlTO^SCALC *f t K I ID 
AC 1 1 Tri J.IMS«>PrAlT0P*2.0l/N.0 
DO IUOO n:l,N«ACH 
IFINPPINT.CQ.OIGO TO 10 
WRITCI6,<tl 

WP1TCI6.SIXNIIKI .PIIKI.PBOPIIKI.SPIIMI 
WBl Hit .71 
10 CONTINUE 

00 100 1 : 1 . NN02 
WRI TC 16 .61 
Jj:0 

NPI I ,k I to 
00 SO JM.JNAI 

1MTHCTLII.J1.LC .O.OIGO TO 50 
ASTAR-AEXIT/AOASI 1. J» 

P ON A X : 1 .BePWNAXPSQRT ITO I /AS TAP 
IFIPONAX.GT.PCNAXIPONAxrPCNAX 
OCNONI I . J I : INC I I .JIppAPGANA JPPB 
CP INC F OL LOWING STATCNCNTS CALCULATE SPNAX 
PCU0P1 1 I .Jl :PONA»/PI IN | 

PCUOPB I I ,J) iPCUOPI I 1 , Jl /f GOP I IN I 
CP NACM NUNBCR ON THC JE T BOUNDARY • NAI1NUN 
CALL NACMIGANAJi PC UOPB I l.JI , XNJNAXII, Jl I 
CP JL T BOUNOARV ISCNTR0P1C C Xp ANSI ON ANGLE - NAXlNUN 
CALL PRNOIGAHAJ, XNJHAXI l.JI .PNJNAXI I, Jl I 
C* PLUNC BOUNOAHY INITIAL EXPANSION ANGLE - NAXlNUN 
OCL JNI I I ,JI:THCUI I . Jl» PNJNAXI I , JI-PNM l.JI 
SPNAX I 1 , Jl: XNJHAXI I. JIPOCLJNXII .JI/CCNONII .Jl 
C* Th[ FOLLOWING STATCNCNTS CALCULATE SPHIN 
PCLOPBI I.JMPCOPII l.JIPPCPBNN 
C* NACH NUNBCR ON THC JC T BOUNDARY - NININUN 
CALL NACHIGANAJ.PCLOPBI l.JI .XNJNINII .Jl I 
C* JC T BOUNOARV ISCNTRQPIC EXPANSION ANGLE - NININUN 
CALL PR NDIGANAJ. XN JNI N I l.JI .PNJNINII, Jl I 
CP PLUNC BOUNOARV INITIAL EXPANSION ANTlC - NININUN 
DELJNNI I,JI:THCTlI 1.JI*PNJN1NII.JI-PNCI1.JI 
SPNINII , JlrxPJNINI I.JIP0CLJNNI1 .JI/OCNONII.JI 

I F I SP1 IK I.GT. SPNINII , Jl ♦3.0.AN0.SP1 IN I .L T .SPNAX 1 1 , J I -5 .0169 TO «0 
GO TO 50 
NO CONTINUE 
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Subroutine SHRNGE (Cont'd) 

IFtNPRINT.CC.OIGO 10 NS 

HPITCIt,)IXMC II, Jl .TMCTLU, JI.SPMINI! ,jl .SPMAXII ,JI,P0MAX 
NS CONTINUE 
jj: JJ*1 

NPI I ,« I SNPI I ,« I • | 

IMC0'I.JJ,MI:XNCII,J> 

TMCTLOt I,JJ,X IsTHETLII, J| 

50 CONI I hUC 
ICO CONTINUE 
1000 CONTINUE 
RETURN 
f NO 


subrcut INC TEST 

DIMENSION Hirer I SI , TMLQF 151 , AOASFIS) ,DST ARE IS I, SAVE II I 
COMMON/ CQ1 /GA «AJ ,GA MAE , GAMA I ,PATM,A ,B 
COMMON/C02/PCMAX ,PEPBMN,PC 

COMMON/ CO*/ SP 11 20 1 , POOP 1 1 20 I , XM 1 1 20 I • PI 1 20 1 
C CM MON/ C 01 /SC AtC «PC 1 1 TO » WMA * * TO 

COMMON/ TEST l/*MC L,XMER, THLO . IHL T ,TMUT , THUS, UMAX 
COMMON/ TEST 2/ SNO If 

1 F ORM A T I 1 MO (IHi* POWER SHEEP OPERATING CHARACTERISTICS FOR EACH OF T 
THE DESIRED Ft IGH 7 CONDITIONS TO BE SIMULATED IN HIND T UNNCt * /2l » *T 
It STS IS PRESENTED BELOH FOR REPRESENTATIVE N022LCS THAT EXIST IN T 
SHE FINAL ENVELOPE LEFINED ABOVE. •» 

2 rORMATI 1H0.I1H ITSUB HNC BE 10. SI 

S roRMAIMHC.lX, • N02 2LE MCXIT LIP ANGLE M I NF PCMIN 

I SPM1N PCNOM SPNOM PC«AX S» MAX HlNO TUNNEL DSTA 

IR Of XIT •/ IOS x, ’PINT • I 

N f ORHAf I 1H0I 

S FORMA* I lX.IS.'X ,9E 10.S, IX , 3t 10.SI 

t FCRMATI 1HC,IX,'MJETF IS SUBSONIC AND TEMPORARILV SET r 1.10*1 
C* THIS SUBROUTINE DETERMINES ThE POHER ShEEP OPERATING ChARACTCRIS- 

C* TICS or REPRESENTATIVE N022LES THAT EXIST IN THE FINAL ENVELOPE 

C* FOR EACH OF THE DESIRED FLIGHT CONDITIONS TO BC SIMULATED IN HIND 

C* TUNNEL TESTS. 

HRI TE lb , 1 I 
HRITEI6.J) 

X "EOF I 1 !:lXMEL*XM£RI/2.0 
X M£ 0 F 12 I :xMEOF I 1 I 
xmeoF I3I:xmeDF 1 1) 
xmeofi* cxmel 

XMEOF IS UXMEP 

TMLDFI2l:lTHLT*THUTI/2.D 

THLDF I 3 I :i THLB* T hus 1/2. 0 

THLDFI1 1 : n hlCF I 2 I ♦ IhlDM * I 1/2.0 

THLDF I* i:iTHLE*THLTI/2.0 

1HL0FI5»;lTHUPrTHUTI/2.0 

CXITD-SCALt*PfKlTO 

AC X I T - I 3. 14 1S9*C Xl TD**2 .01/4.0 

DO 1001 l:i,NN02F 

HRITEI6.4I 

C* NC22LI EXIT PLANE AREA RATIO AC/A* 

CALL AUATIOIGAHAE , XMEOF III, AOASF II I I 

ASTAR-AEXIT/AOASF I 1 1 

DSTARFI I I :£ XI TO /SORT I AO ASF I III 

C* MAXIMUM ALLOHABLE CHA1BCR PRESSURE FOR A GIVEN N022LE SUBJECT TO 

C* THE TEST FACILITY MASS FLOH CAPABILITY 
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I 


Subroutine TEST (Cont'd) 

POMAXri .aBPNMAXPSORT ITO l/ASTAR 

oenohfixmeof i iippapgama jppb 
THE FOLLOWING STATEMENTS CALCULATE SPMIN 
N022LE LIP ISTNTROPIC EXPANSION ANGLE 
CALL PRNOIGAMAE, XMEOF II I.PMEF I 
N022LE EMIT PLANE PRESSURE RATIO fC/P C 

CALL PR A T 10 1 GAMA E , XMEOF 1 1 I « PCOPEF I 
PCLPBFrPCOPEF pPEPBMN 

MACH NUMBER ON JET BOUNDARY - MINIMUM 
CALL MACHIGAMAJ'PCLPBF'KHJMNFi 

JET BOUNOAR V ISENTROPIC EXPANSION ANGLE - MINIMUM 
CALL PRNOIGAMAJ, XMJMNF,PMJMNFI 

PLUMC BOUNDARY INITIAL EXPANSION ANGLE - MINIMUM 
DC JMNF: THLOF III «PM JMNF-PMEF 
SPM I NF I XM JMNF POE JMNF /DC NOME 
DO 1000 k:1,kmAX 

THE FOLLOWING STATEMENTS CALCUI ATE SPMAX 
PCOPIF=POMAX/PIIKI 
PCOPBFiPCOPIF /PBOPIINI 

MACH NUMBER ON JET BOUNOARY - MAXIMUM 
CALL MACMIGAMAJ.PCOPBF, XMJMXFi 

JET BOUNDARY ISENTROPIC EXPANSION ANGLE - MAXIMUM 
CALL PRNOIGAMAJ, XMJMXF t PMJMXF I 

PLUMC BOUNDARY INITIAL EXPANSION ANGLE - MAXIMUM 
OE JMXFrTHLDE l J MPMJMXF - PMEF 
SPMAXFSXMJMXF POE JMXF/OENOMF 

MINIMUM CHAMBER PRESSURE REQUIRED TO PREVENT FLOW SEPARATION 
POMIN:PCLPBF PPBOPI IX I PP I IN I 

THE FOLLOWING STATEMENTS CALCULATE THE RCGUIRCD CHAMBER PRESSURE 
10 SIMULATE THE DCSIREO FLIGHT CONDITIONS 
XMJF : 1 . 1PXMC0FI I I 
ifcO CONTINUE 

save m = i.o 

SAVEl?l:.10 
T0LS.P5 
170 CONTINUE 

CP JET BOUNOARY ISENTROPIC EXPANSION ANGLE 
CALL PRNOIGAMAJ, XMjF ,PM JF » 

C* PLUME BOUNDARY INITIAL EXPANSION ANGLE 

DEL JF:| SPI IN I /XHJP IPOENOMF 
CP NO? 2 LE LIP ANCLE 

THE TLP:DEIJF ♦PMEF-PMJF 
DTHETL:ThLDFI I l-THETLP 
CALL I T SUB ID THE TL ,XMjF,SAVC ,TOL ,1991 
I I = I F I X I SAVE I 1 ) • . 1 0 I 
N 1 1 1 F I X I SAVE m».IOI 
IF (XMJF .GT.l .LIGO TO II 
XMJF :| . 10 
WR1 TE 16 ,6 I 
GO TO 180 
11 CONTINUE 

GO TCI170,17C,17r,17C,18C,17SI,II 
175 CONTINUE 

W R1 T 1 1 6 , 2 I TKL DF I I » , THE TLP, XM JF .DEL JF 
160 CONTINUE 

CP PLUME BOUNDARY PRESSURE RATIO PC/PB 
CALL PRATIO IGAMA J,XMJF ,PCOPBF I 
PONOMrPCOPBF APBOPIINIpPIINI 

WRITCI6,SII,XMC0FIII,THL0FI II , XMI IK I , POMIN , SPMINF ,PONOH ,SP1 IN I , 
«POMAX,SPMAXF,PIIKl,OSTAPFII I, EX I TO 
1 COD CONTINUE 
1001 CONTINUE 
RETURN 
END 
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CP 

CP 

CP 

CP 


CP 

CP 
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CP 
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SUBftOUT INC IVLIN (SHIN, SMXX ,UU , VB , V T , OV , MM , 00 ,NV t 

COMMON XN,VN,2N 

B£*L INC 

INTEGER UUtOO 

NNSO 

1FIUU.EQ.5IG0T065 

ir (uu.eq.ii go to to 

VT:0.0 

ve:o .n 

XINC:(SHAX-SMJM /10.0 

1MXINC.LT. 1.0. AND. XINC.GE. 0.061 GO TO 30 
IMIINC.LT. 0.06) GO TO 65 

00 10 1=1. SO 
NN- I - 1 

IN- . 5 

IF(2.*10.**( 1*1 I .GT.XINCI GO TO 20 
IN: 1 . 

IMS. *10. **11-11. GT.XINCI GO TO 20 
X N: 2 • 

IM 1 .*1 0.** I I I.G T.XINC I GO TO 20 
1L C0NT1NUC 
2 C INC:XN* 10.**NN 
GO TO 60 

jc oiff:xinc*io.d 

Nn:NN« 1 
INC:.P? 

irio irr . gt. i .si inc:.os 

IMOirr .GT. 3.01 1NC= . 1 
IMOirr .GT.7.PI INC:. 2 
tc continue 

YT:SNXX . AMOD ( *BS I SHXXI , (INC *5. II 

V !* :SMIN ~AMOD I IBS <SMINI,(INC*5.II 

IMSHAX.GT.O.OI YT:SMAX«xeS UMOOiXBSI SMXX I , I INC*5 . I I - ( INC*5 . 1 1 
IMSMIN.LT. 0.01 VB :SM I N -AB S ( XMOO ( »BS ( SM IN I , ( INC* S . I I - ( INC*S . I I 
OV: I NC 
mm: S 
00=10 
NV:NN« I 
GO TO 100 
65 CONTINUE 
X:$M«X* 10. 
x:ir 1*( XI 
VT:A/10.*0. 1 
X :SMIN* 10. 

X : I r I * ( X I 

V [• : X / 1 0 . 

1 NC : 0 . 0 1 
GOT 0 60 

TG DO 80 1:1.20 
C LOG SCALES 

C the SMXLLCST VXLUE FOR SMXX IS 10**-4 

NNN: I -5 

80 IF (SMXX ,LT . 10.**NNNIG0 TO 9 
9 V T: 1 Q.0**NNN 
00 90 1=1,20 
NNN: 10- I 

C THE LARGEST VALUE FOR SMIN IS 10**9 

90 IF(SMIN.GT.IO.O**NNNIGO TO 91 

91 VB :i 0 .0 **NNN 
0 V : 1 . r 

MM: 1 
00: 1 
NV--2 

call smxvv (o.n 

ICC CONTINUE 
RETURN 
END 
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